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INTRODUCTION 

Triple  negative  breast  cancers  (TNBCs)  account  for  ~10-15%  of  breast  cancers  and  are 
associated  with  poor  clinical  outcomes  and  worse  disease-free  and  overall  survival  [1],  Thus, 
there  is  a  critical  need  to  identify  new  therapeutic  targets  for  this  breast  cancer  subtype.  TNBCs 
lack  of  expression  of  three  receptors  that  are  commonly  used  to  determine  the  prognosis  and 
treatment  options  for  patients:  estrogen  receptor  (ER)  a,  progesterone  receptor,  and  human 
epidermal  growth  factor  receptor  2  (HER2).  A  second  estrogen  receptor  (ER(3)  has  been  detected 
in  some  TNBCs  [2-4]  and  exogenous  expression  of  ERp  inhibits  the  growth  of  breast  cancer 
cells  [5-8].  Given  the  anti-proliferative  activity  of  ERp,  this  receptor  may  be  an  effective 
therapeutic  target  for  triple  negative  breast  cancers  that  express  ERp.  The  hypothesis  underlying 
this  project  is  that  ERfi  is  an  effective  therapeutic  target  and  ERf  selective  ligands  can  inhibit 
the  growth  of  triple  negative  breast  cancers  when  ER^  is  expressed. 

Three  aims  were  originally  proposed  to  test  this  hypothesis:  1)  characterize  the  growth 
inhibotry  effects  of  ERP-selective  ligands  in  TNBC  cells;  2)  identify  new  ERP-selective  ligands 
using  high  throughput  screening;  and  3)  determine  the  effects  of  ERP-selective  ligands  on  the 
growth  of  TNBC  in  xenograft  models.  Based  on  the  data  collected  and  feedback  from  both  the 
committee  and  the  reviewers  of  the  original  proposal,  some  of  the  original  aims  and  experiments 
were  re-designed  and  a  revised  statement  of  work  (SOW)  was  approved  in  2012.  This  revised 
SOW  incorporated  2  new  aims:  1)  characterize  the  growth  inhibitory  effects  of  ERp  and  identify 
ERp  target  genes  using  RNA  sequencing  (RNA-seq);  2)  assess  ERp  expression  in  a  cohort  of 
TNBC  clinical  samples.  This  report  will  detail  the  results  of  each  aim  from  the  original  SOW  and 
the  revised  SOW.  Key  research  accomplishments  and  reportable  outcomes  will  also  be 
described. 

BODY 

Estrogen  signaling  is  primarily  mediated  by  two  estrogen  receptors  (ERs),  ERa  and  ERp. 
Despite  sharing  several  structural  similarities,  ERa  and  ERp  have  many  unique  features  with 
respect  to  ligand  binding  and  gene  regulation.  ERa  and  ERp  are  expressed  in  a  variety  of  tissues 
and  have  both  overlapping  and  distinct  tissue  distributions,  and  there  are  several  ERp  isoforms 
that  can  modulate  the  estrogen  response.  ERs  can  regulate  gene  expression  in  both  ligand 
dependent  and  independent  manners.  ER  ligands  are  structurally  diverse,  and  some  ligands  show 
selectivity  for  ERa  or  ERp  due  to  differences  in  the  ligand  binding  pockets  of  the  two  receptors. 
Ultimately,  ligands  can  elicit  tissue-  and  cell  type-specific  responses,  as  exemplified  by  selective 
estrogen  receptor  modulators  (SERMs)  such  as  tamoxifen,  which  is  used  to  target  ERa  for  breast 
cancer  treatment.  ERa  is  expressed  in  -70%  of  breast  cancers,  but  -10-15%  of  breast  cancers  are 
classified  as  triple  negative  breast  cancers  (TNBCs).  TNBCs  lack  therapeutic  targets  and  there  is 
a  need  to  identify  new  treatment  strategies  for  this  particularly  aggressive  breast  cancer  subtype. 
ERp  is  expressed  in  a  subset  of  TNBCs,  and  previous  in  vitro  data  suggested  that  ERp  is  growth 
inhibitory  in  ERa-negative  breast  cancers.  Based  on  the  results  of  previous  studies,  it  was 
hypothesized  that  ERp  may  be  a  therapeutic  target  in  a  subset  of  ERa-negative  breast  cancers, 
including  some  TNBCs. 
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Aim  1:  Characterize  the  inhibitory  effects  of  ERp-selective  compounds  in  TNBC 

The  goals  of  this  aim  were  to:  1)  determine  the  effeets  of  eosmosiin  and  liquiritigenin  on 
the  growth  of  Hs578T-ERp  eells;  2)  identify  additional  models  of  ERP-positive  TNBC;  and  3) 
identify  downstream  targets  of  ERP  that  mediate  growth  inhibitory  effects.  Much  of  the  work 
towards  this  aim  has  been  published  [9,  10]  (Appendix  B)  and  was  described  in  the  previous 
annual  reports.  Briefly,  eosmosiin  and  liquiritigenin  inhibited  the  growth  of  Hs578T-ERp  cells  (a 
TNBC  cell  line  with  inducible  ERP  expression)  only  when  ERP  was  expressed.  This  inhibitory 
effect  is  likely  due  to  a  cell  cycle  arrest  because  these  compounds  did  not  induce  apoptosis  (data 
not  shown)  and  17b-estradiol  (E2),  the  most  prevalent  estrogen  in  women  of  reproductive  age, 
was  shown  to  induce  a  cell  cycle  arrest  in  Hs578T-ERp  cells  [11]. 

Next,  these  results  were  confirmed  in  another  model  of  ERP-positive  TNBC.  Because 
additional  models  of  TNBC  with  detectable  endogenous  ERp  could  not  be  identified,  a 
tumorigenic  TNBC  cell  line  (MDA-468)  was  engineered  to  express  ERP  after  doxycycline  (Dox) 
treatment.  These  MDA468-ERP  cells  were  described  in  the  last  annual  report,  and  since  then 
have  been  described  extensively  in  a  recent  publication  [10]  (Appendix  B).  Briefly,  E2  induced  a 
GO/Gl  cell  cycle  arrest  in  MDA468-ERP  cells  only  when  ERP  was  expressed,  and  surprisingly, 
ERp  expression  alone  induced  cell  cycle  arrest,  although  to  a  lesser  degree.  These  results 
suggested  that  these  cells  could  be  used  to  identify  ligand-dependent  and  ligand-independent 
ERp  target  genes  in  TNBC  cells. 

In  order  to  identify  downstream  targets  of  ERp  that  mediate  growth  inhibition,  RNA-seq 
was  performed  in  MDA468-ERP  cells  in  the  presence  and  absence  of  E2  [10].  This  approach 
allowed  the  first  global  identification  of  ligand-dependent  and  ligand- independent  ERp  target 
genes  in  the  absence  of  ERa,  thereby  providing  a  rich  resource  for  the  scientific  community. 
Since  this  work  is  published  [10],  only  a  few  of  the  results  will  be  highlighted  in  this  report. 

Eirst,  ERp  induced  the  expression  of  the  cyclin  dependent  kinase  inhibitor  p21,  which  mediates 
the  progression  of  the  cell  cycle.  Second,  several  common  ERp  target  genes  were  identified  in 
Hs578T-ERp  and  MDA468-ERP  suggesting  that  the  receptor  can  mediate  a  common  set  of  genes 
in  different  TNBC  cell  lines.  Einally,  in  collaboration  with  a  biostatistician,  ERp  expression  was 
found  to  be  associated  with  the  expression  of  several  of  its  target  genes  in  TNBC  gene 
expression  data  from  the  Cancer  Genome  Atlas  database.  Overall,  this  is  the  most  comprehensive 
study  to  date  regarding  target  gene  regulation  by  ERp.  In  addition,  this  study  allowed  me  to  build 
a  bioinformatics  collaboration  and  develop  the  skills  necessary  to  perform  large-scale  genomics 
studies.  These  skills  will  be  critical  for  my  future  career  in  breast  cancer  research. 

Original  Aim  2:  Identify  new  ERp-selective  ligands  using  high  throughput  screening  (HTS) 

The  goals  of  this  aim  were  to:  1)  optimize  luciferase  assays  with  Hs578T-ERpEuc,  2) 
screen  additional  Hs578T-ERaEuc  cells,  and  3)  conduct  a  screen  and  counterscreen  with  small 
molecule  libraries  using  Hs578T-ERaLuc  and  Hs578T-ERpEuc  cells.  Although  significant  effort 
went  into  the  high  throughput  screening  optimization,  this  aim  was  ultimately  deleted  because  of 
the  variability  in  the  screening  assay.  As  described  in  the  previous  annual  reports,  the  assay 
could  not  be  optimized  for  a  384  well.  However,  these  reporter  cell  lines  were  highly  sensitive  to 
estrogenic  ligands,  and  this  part  of  the  project  ultimately  resulted  in  two  first- author  publications 
[9,  12]. 
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Original  Aim  3:  Determine  the  effects  of  ERp-selective  ligands  on  the  growth  of  TNBC  in 
xenograft  mouse  models 

The  original  goals  of  this  aim  were  to:  1)  graft  Hs578T-ERaLuo  and  Hs578T-ERpLuc 
cells  into  nude  mice;  2)  inject  mice  with  ERP  ligands  and  monitor  tumor  growth  using 
bioluminescent  imaging;  and  3)  assess  luciferase  expression  and  tumor  histology  after 
treatments.  Although  several  attempts  were  made  to  graft  Hs578T-ERaEuc  and  Hs578T-ERpEuc 
into  nude  mice,  palpable  tumors  did  not  consistently  form  and  this  part  of  the  aim  was  revised. 

As  an  alternative  approach,  MDA468-ERP  cells  were  labeled  with  luciferase  and  injected 
into  the  fatpad  of  nude  mice  to  track  the  growth  of  the  tumor  after  treatment  with  Dox  and/or  E2. 
The  results  of  this  study  are  published  [10],  and  the  growth  inhibitory  effects  of  ERp  were 
confirmed  in  vivo. 

Revised  Aim  3:  Determine  the  effects  of  ERp  expression  and  activation  on  the  growth  of 
TNBC  in  xenograft  mouse  models  and  assess  ERp  expression  in  TNBC  clinical  samples 

The  first  part  of  this  aim  is  described  above.  In  order  to  address  concerns  of  the  reviewers 
and  thesis  committee,  a  second  part  of  the  aim  was  incorporated  into  the  revised  statement  of 
work.  A  major  weakness  of  the  original  proposal  was  that  it  is  not  widely  accepted  that  ERp  is 
expressed  in  TNBCs.  Many  studies  aimed  at  assessing  ERp  expression  have  yielded  inconsistent 
results  (reviewed  in  [13]).  In  order  to  address  the  issues  of  antibody  specificity, 
immunohistochemistry  (IHC)  protocols  for  ERp  detection  were  optimized  using  MDA468-ERP 
xenograft  tissues  in  which  the  expression  of  ERp  was  regulated  by  Dox  exposure.  As  shown  in 
Eigure  1,  ERp  was  specifically  detected  in  the  xenograft  tissues  obtained  from  mice  exposed  to 
Dox  (+ERP)  and  this  signal  was  eliminated  by  preabsorbing  the  antibody  with  peptide 
corresponding  to  ERp.  These  IHC  protocols  were  then  applied  to  two  cohorts  of  TNBC  from  the 
Marshfield  Clinic  in  collaboration  with  the  Translational  Research  Initiatives  in  Pathology 
(TRIP)  lab  and  Dr.  Kari  Wisinski  (manuscript  in  preparation)  (Eigure  2).  In  order  to  objectively 
score  ERp  expression  in  the  samples,  the  VECTRA^'^  multispectral  imaging  instrument  was 
used  to  quantify  the  ERp  signal  in  the  nuclear  and  cytoplasmic  compartments.  As  shown  in 
Eigure  3,  there  was  a  close  relationship  between  the  percent  positive  nuclei  and  the  mean  optical 
density  (OD)  for  ERp  detection.  There  was  also  a  close  relationship  between  the  cytoplasmic  and 
nuclear  ERp  expression  (Eigure  3B  and  3C). 

The  tissues  were  also  stained  for  Ki67,  an  indicator  of  proliferation,  to  assess  whether 
ERp  expression  was  associated  with  proliferation.  Surprisingly,  there  was  a  significant  positive 
correlation  between  ERp  and  Ki67  (Eigure  4 A  and  4B).  A  survival  analysis  revealed  that  there 
was  no  association  between  ERp  expression  and  survival  in  this  cohort  of  TNBCs  (Eigure  5  A 
and  5B).  These  results  are  currently  being  prepared  for  publication  and  will  surely  be  published 
by  the  end  of  the  year  because  this  is  the  first  objective,  quantitative  analysis  of  ERp  expression 
in  TNBC.  Although  there  was  a  discrepancy  between  the  antiproliferative  effects  of  ERp  in  vitro 
and  the  positive  association  between  ERp  and  Ki67  in  clinical  samples,  this  part  of  the  project 
provides  a  foundation  on  which  ERp  IHC  should  be  performed  in  future  studies  with  larger 
cohorts.  In  addition,  this  aspect  of  the  project  provided  experience  in  bridging  basic  and  clinical 
research  that  is  necessary  for  successful  translational  research.  I  learned  the  challenges 
associated  with  collaborating  and  communicating  with  clinical,  statistical  and  basic  researchers, 
and  this  invaluable  experience  will  help  me  successfully  drive  these  translational  projects  in  the 
future. 
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KEY  RESEARCH  ACCOMPLISHMENTS: 

•  Characterization  of  growth  inhibitory  effects  of  cosmosiin  and  liquiritigenin  in  Hs578T- 
ERp  cells 

•  Generation  of  MDA468-ERP  inducible  cell  lines 

•  Characterization  of  the  sensitivity  and  selectivity  of  Hs578T-ERpEuc  and  Hs578T- 
ERaEuc  cells 

•  Optimization  of  HTS  assay  for  Hs578T-ERpEuc  in  96  well  format 

•  Determined  that  ERP  inhibits  MDA468-  ERP  cell  growth  by  inducing  cell  cycle  arrest 

•  Demonstrated  that  ERp  regulates  p21,  a  key  regulator  of  G1  arrest 

•  Developed  a  retroviral  delivery  system  to  express  ERp  in  additional  TNBC  cell  lines 

•  Identification  of  ERp  targets  using  RNA  sequencing  using  RSEM  and  DESeq 

•  Identified  several  functions  and  pathways  that  are  enriched  in  ERb  target  genes  and 

•  validated  the  ERp  target  genes  using  quantitative  PCR 

•  Developed  MDA468-  ERp  -EEuc  cells  for  xenografts  experiments  and  in  vivo  imaging 

•  Assessed  the  expression  of  ERp  in  a  cohort  of  ERa-negative  breast  cancer  samples 

•  Determined  the  associations  between  ERp  expression  and  clinicopathologic 
characteristics  in  a  cohort  of  TNBC 

REPORTABLE  OUTCOMES: 

Conferences  Attended: 

Society  of  Toxicology  Annual  Meeting,  2012 

Gordon  Research  Conference  on  Hormone  Action  in  Development  and  Cancer,  2011 

Awards: 

RSESS  Society  of  Toxicology  2012  Annual  Meeting  Travel  Award,  2012 

Publications: 

Shanle  E,  Onitilo  A,  Huang  W,  Kim  K,  Zang  C,  Engel  J,  Xu  W,  Wisinski  K.  Prognostic 

significance  of  full  length  estrogen  receptor  beta  expression  in  Stage  I-III  triple  negative 
breast  cancer  .  In  preparation. 

Yarger  JG,  Babine  RE,  Bittner  M,  Shanle  E,  Xu  W,  Hershberger  P,  Nye  SH,  2012.  Structurally 
similar  estradiol  analogs  uniquely  alter  the  regulation  of  intracellular  signaling  pathways. 
J  Mol  Endocrinol.  50,  43-57. 

Sievers  C*,  Shanle  EK*,  Bradfield  C,  Xu  W,  2012.  Differential  action  of  monohydroxylated 
polycyclic  aromatic  hydrocarbons  with  estrogen  receptors  a  and  p.  Toxicol  Sci  [Epub 
ahead  of  print].  *  Authors  contributed  equally  to  this  work 

Powell  E,  Shanle  E,  Brinkman  A,  Ei  J,  Keles  S,  Wisinski  K,  Huang  W,  Xu  W.  Identification  of 
Estrogen  Receptor  dimer  selective  ligands  reveals  growth-inhibitory  effects  on  cells  that 
co-express  ERa  and  ERp.  PLoS  ONE,  e30993. 

Shanle  E,  Xu  W.  Eunction,  expression,  and  detection  of  estrogen  receptor  isoforms  in  normal 
and  malignant  tissues,  in:  Chen,  G.  (Ed.),  Estrogen  Receptors:  Mechanisms,  Structure 
and  Role  in  Disease.  Nova  Science  Publishers  Inc.,  New  York,  ISBN  978-1-62257-180-2. 
(Book  Chapter) 

Shanle  E,  Xu  W.  Generation  of  stable  reporter  breast  cancer  cell  lines  for  the  identification  of 
ER  subtype  selective  ligands.  Biochem  Pharmacol  82,  1940  (2011). 
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Shanle  E,  Xu  W.  Endocrine  disrupting  chemicals  targeting  estrogen  receptor  signaling: 
Identification  and  mechanisms  of  action.  Chem  Res  Toxicol  24,  6  (2010).  (review) 


Poster  Presentations: 

Blanke,  K.,  Shanle,  E.  Exploring  Toxicology:  Designing  learning  goals  and  evaluation  strategies 
for  outreach  activities;  Accepted  for  oral  presentation  in  the  K-12  Toxicology  Outreach 
Activities:  Regional  Chapter  Successes  and  Resources  Informal  Specialty  Section.  Society 
of  Toxicology  Annual  Meeting,  March  13,  2013. 

Shanle,  E,,  Xu,  W.  Stable  reporter  cell  lines  for  the  identification  of  subtype  selective  estrogenic 
ligands.  Society  of  Toxicology  Annual  Meeting,  March  11-15,  2012.  San  Erancisco,  CA. 
Poster  Presentation. 

Shanle,  E,,  Xu,  W.  Stable  reporter  cell  lines  for  the  identification  of  subtype  selective  estrogenic 
ligands.  Gordon  Conference  for  Hormone  Action  in  Development  and  Cancer,  July  30- 
August  4,  201 1.  Bryant  University,  Smithfield,  RI. 

Shanle,  E.,  Powell,  E.,  Eastlund,  V.,  Xu,  W.  Growth  inhibitory  effects  of  natural  phytoestrogens 
on  breast  cancer  cells  mediated  by  ERp.  Great  Eakes  Nuclear  Receptor  Conference,  October 
22-23,  2010.  University  of  Michigan,  Ann  Arbor,  MI. 

CONCLUSION 

At  the  completion  of  the  Department  of  Defense  Breast  Cancer  Research  Program  Pre- 
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Appendix  A:  Figures  and  Tables 

Figure  1:  Xenograft  tissues  with  inducible  ERP  expression  are  useful  for  optimizing  ERP  IHC. 
MDA468-ERP  cells  were  injected  into  the  mammary  fat  pads  of  nude  mice.  After  tumors 
formed,  mice  were  treated  with  either  vehicle  (1%  sucrose)  (A,  B,  C)  or  Dox  (D,  E,  E).  IHC  was 
performed  with  the  PAl-313  ERp  antibody  (A,  D).  Eor  controls,  the  antibody  was  pre-absorbed 
with  ERp  peptide  (B,  E)  or  the  primary  antibody  was  excluded  entirely  (C,  E).  The  brown 
staining  that  indicates  reactivity  toward  ERp  is  only  observed  in  tissues  from  mice  exposed  to 
Dox  (+ERP)  (A).  Scale  bars  =100  pm. 
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Figure  2:  Immunohistochemistry  staining  for  ERP  in  TNBCs  from  the  Marshfield  Clinic.  A  and 
B)  Representative  images  of  ERP-negative  tumor  tissues  (0.33%  and  2.6%  1+  or  higher  nuclei, 
respectively).  C  and  D)  Representative  images  of  ERP-positive  tumor  tissues  (66.2%  and  73.2% 
1+  or  higher  nuclei,  respectively).  Scale  bars  =  200  pm. 
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Figure  3:  Relationships  between  ERP  subcellular  localization  and  scoring  strategies  in  the 
Marshfield  cohort.  A)  A  comparison  of  the  nuclear  mean  OD  and  percent  ERP-positive  nuclei. 
B)  Relationship  between  the  nuclear  and  cytoplasmic  ERP-positive  cells.  C)  Relationship 
between  the  nuclear  and  cytoplasmic  mean  OD  values. 
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Figure  4:  ERP  expression  is  correlated  with  Ki67  expression.  A  scatter  plot  showing  the 
relationship  between  the  percentage  of  Ki67-positive  nuclei  (1+  or  greater)  and  the  mean  OD  for 
(A)  nuclear  ERP  expression  or  (B)  cytoplasmic  ERP  expression  as  detected  using  the  PAl-313 
antibody.  There  was  a  moderately  significant  positive  correlation  between  the  two  variables 
(Spearman  correlation,  p<0.0001  for  both  comparisons). 
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Figure  5:  ERP  expression  is  not  signifieantly  assoeiated  with  disease-free  survival  or  overall 
survival  in  the  Marshfield  eohort.  Kaplan-Meier  graphs  for  disease-free  survival  (A)  and  overall 
survival  (B)  with  respeet  to  ERP  low  (<20%  nuelear  staining,  n  =  42)  and  ERP  high  (>20% 
nuclear  staining,  n  =  28)  expression. 
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Abstract 

Estrogen  signaling  is  mediated  by  two  estrogen  receptors  (ERs),  ERa  and  ERp,  which  have  unique 
roles  in  the  regulation  of  breast  cancer  cell  proliferation.  ERa  induces  proliferation  in  response  to 
estrogen  and  ERp  inhibits  proliferation  in  breast  cancer  cells,  suggesting  that  ERp  selective 
ligands  may  be  beneficial  for  promoting  the  anti-proliferative  action  of  ERp.  Subtype  selective 
ligands  can  be  identified  using  transcriptional  assays,  but  cell  lines  in  which  ERa  or  ERp  are 
independently  expressed  are  required.  Of  the  available  reporter  cell  lines,  none  have  been 
generated  in  breast  cancer  cells  to  identify  subtype  selective  ligands.  Here  we  describe  the 
generation  of  two  isogenic  breast  cancer  cell  lines,  Hs578T-ERaLuc  and  Hs578T-ERpLuc,  with 
stable  integration  of  an  estrogen  responsive  luciferase  reporter  gene.  Hs578T-ERaLuc  and 
Hs578T-ERpLuc  cell  lines  are  highly  sensitive  to  estrogenic  chemicals  and  ER  subtype  selective 
ligands,  providing  a  tool  to  characterize  the  transcriptional  potency  and  subtype  selectivity  of 
estrogenic  ligands  in  the  context  of  breast  cancer  cells.  In  addition  to  measuring  reporter  activity, 
ERp  target  gene  expression  and  growth  inhibitory  effects  of  ERp  selective  ligands  can  be 
determined  as  biological  endpoints.  The  finding  that  activation  of  ERp  by  estrogen  or  ERp 
selective  natural  phytoestrogens  inhibits  the  growth  of  Hs578T-ERp  cells  implies  therapeutic 
potential  for  ERp  selective  ligands  in  breast  cancer  cells  that  express  ERp. 


Keywords 

Estrogen  receptors;  subtype  selectivity;  phytoestrogens;  breast  cancer 


1.  Introduction 

Estrogens  regulate  mammary  gland  growth  and  differentiation,  ovary  and  uterus  maturation, 
and  bone  homeostasis  [1].  The  physiological  effects  of  estrogens  are  primarily  mediated  by 
two  estrogen  receptors  (ERs),  ERa  and  ERp.  Because  of  the  broad  range  of  ER  target  tissues 
and  the  ligand  dependent  activity  of  the  receptors,  synthetic  and  natural  estrogens  hold 
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therapeutic  promise  in  selectively  targeting  ERs.  Therapies  aimed  at  preventing  ERa 
transcriptional  activation  are  currently  used  for  breast  cancer  treatment  and  osteoporosis 
prevention  [2].  Though  ERp  is  not  currently  a  therapeutic  target,  accumulating  evidence 
suggests  an  anti-proliferative  role  for  ERp  in  breast  cancer  [3].  In  the  mammary  gland,  ERa 
and  ERp  play  opposing  roles  in  regulating  growth  and  differentiation  in  response  to 
estrogens;  ERa  promotes  proliferation  while  ERp  inhibits  ERa-mediated  proliferation  [4-6]. 
Because  the  anti-proliferative  action  of  ERp  may  be  enhanced  by  ligand-dependent 
activation,  the  paradigm  of  ER  targeted  therapies  is  expanding  towards  the  development  of 
ER  subtype  selective  ligands  [7]. 

Though  ERa  and  ERp  share  many  structural  and  transcriptional  features,  ligands  can  display 
subtype  selectivity.  In  classical  ligand  dependent  transcriptional  activation,  the  receptors 
dimerize  upon  ligand  binding  and  undergo  conformational  changes  to  allow  cofactor 
recruitment.  The  receptors  directly  bind  DNA  most  often  at  estrogen  response  elements 
(EREs),  consisting  of  a  consensus  GGTCAnnnTGACC  sequence.  ERa  and  ERp  have  97% 
identity  within  the  DNA  binding  domains,  and  the  receptors  bind  similar  DNA  sequences 
with  high  affinity.  Genome  wide  binding  studies  in  MCF7  breast  cancer  cells  expressing 
ERa  or  ERp  independently  have  shown  that  ERa  and  ERp  bind  similar  sites  in  response  to 
17P-estradiol  (E2);  -60%  of  ER  binding  sites  contain  full  EREs  and  -25%  contain  half 
EREs  [8]. 

The  ligand  binding  pockets  of  ERa  and  ERp  are  relatively  large,  and  the  receptors  bind  a 
wide  array  of  chemicals.  The  ligand  binding  domains  of  ERa  and  ERp  have  59%  identity, 
and  the  receptors  bind  E2  with  similar  affinities.  Despite  similarities  in  their  ligand  binding 
domains,  several  ligands  have  modest  selectivity  for  ERa  or  ERp  [9],  and  some  synthetic 
ligands  maintain  high  selectivity.  For  example,  propyl  pyrazole  triol  (PPT)  is  an  ERa 
selective  agonist  that  displays  a  400-fold  higher  binding  affinity  for  ERa  compared  to  ERp 
[10].  Estrogenic  chemicals  produced  in  plants,  known  as  phytoestrogens,  often  display 
subtype  selectivity  for  ERp.  For  example,  liquiritigenin  is  a  flavanone  derived  from 
Glycyrrhizae  uralensis  that  has  been  shown  to  have  20-fold  higher  binding  affinity  for  ERp 
and  even  greater  selectivity  in  transcriptional  assays  [11].  Compounds  such  as  liquiritigenin 
often  show  low  binding  affinities  relative  to  E2,  and  ERp  selective  ligands  with  higher 
affinity  and  greater  selectivity  are  needed  to  fully  elucidate  the  anti-proliferative  role  of  ERp 
in  breast  cancer. 

Mammalian  cell  lines  have  been  developed  to  enable  screening  for  subtype  selective  ligands. 
FleLa  cervical  carcinoma  cells  have  been  used  to  create  FIELN-ERa  and  FIELN-ERp,  two 
cell  lines  in  which  ERa  or  ERp,  respectively,  are  constitutively  expressed  with  stable 
integration  of  a  luciferase  reporter  downstream  of  an  ERE  [12].  Fluman  embryonic  kidney 
cells,  HEK293,  have  also  been  created  using  a  similar  strategy  in  which  ERa  or  ERp  are 
constitutively  expressed  and  human  placental  alkaline  phosphatase  downstream  of  the 
vitellogenin  ERE  is  stably  integrated  [13].  The  only  available  breast  cancer  reporter  cell  line 
is  T47D-KBLuc  in  which  three  tandem  EREs  upstream  of  a  luciferase  reporter  have  been 
stably  integrated  [14].  Flowever,  identification  of  subtype  selective  ligands  is  prohibited 
because  T47D  cells  express  both  ERa  and  ERp. 

Flere,  we  describe  the  generation  of  two  isogenic  reporter  cell  lines,  FIs578T-ERaLuc  and 
FIs578T-ERpLuc,  that  provide  a  tool  to  characterize  the  transcriptional  potencies  and 
subtype  selectivity  of  estrogenic  compounds  in  the  context  of  breast  cancer  cells.  These  cell 
lines  are  highly  sensitive  to  estrogenic  ligands  and  subtype  selective  ligands  and  can  be  used 
to  validate  ER  transcriptional  activation  by  analysis  of  endpoints  such  as  endogenous  target 
gene  regulation.  Further,  ERp  selective  ligands  are  shown  to  induce  ERP-mediated  reporter 
gene  expression,  endogenous  gene  regulation,  and  growth  inhibition,  suggesting  that 


Biochem  Pharmacol.  Author  manuscript;  available  in  PMC  2012  December  15. 


NIH-PA  Author  Manuscript  NIH-PA  Author  Manuscript  NIH-PA  Author  Manuscript 


Shanle  et  al. 


Page  3 


Hs578T-ERpLuc  cells  may  be  used  to  isolate  ERp  selective  ligands  with  desired  biological 
effects. 

2.  Materials  and  Methods 

2.1  Cell  lines  and  reagents 

Cosmosiin  (apigenin  7-glucoside),  dimethyl  sulfoxide  (DMSO),  E2,  and  diethylstilbestrol 
(DES)  were  obtained  from  Sigma  (St.  Louis,  MO);  DPN,  PPT,  and  ICI  182,780  were 
obtained  from  Tocris  (Ellinsville,  MO);  liquiritigenin  was  obtained  from  Chromadex  (Irvine, 
CA).  Doxycycline  (Dox)  was  obtained  from  Clontech.  Hygromycin  B,  blasticidin  S,  zeocin, 
NaCl,  sodium  dodecyl  sulfate  (SDS),  and  dithiothreitol  (DTT)  were  obtained  from  Research 
Products  International  (Mount  Prospect,  IE).  Triton  X-100  was  obtained  from  Fisher  (Fair 
Lawn,  NJ);  protease  inhibitors  were  obtained  from  Roche  Scientific  (Basel,  Switzerland); 
benzonase  was  obtained  from  Novagen  (San  Diego,  CA).  All  other  chemicals  were  obtained 
from  Sigma  (St.  Louis,  MO). 

Cell  culture  media  were  obtained  from  Invitrogen  (Carlsbad,  CA).  MCF7  and  HEK293  cells 
were  cultured  in  DMEM  +  10%  fetal  bovine  serum  (FBS;  Gemini  Bio  Products,  West 
Sacramento,  CA)  at  37  °C  and  5%  CO2.  Hs578T-ERa  and  Hs578T-ERp  were  previously 
created  by  Secreto  and  coworkers  [15].  These  cells  were  cultured  at  37  °C  and  5%  CO2  in 
DMEM/F12  supplemented  with  L-glutamine,  10%  Tet-system  approved  FBS  (Clontech 
Mountain  View,  CA),  500  mg/L  Zeocin  and  5  mg/L  Blasticidin  S. 

2.2  Generation  of  Hs578T-ERaLuc  and  Hs578T-ERpLuc  reporter  cell  lines 

Stable  reporter  cell  lines  were  created  using  a  modified  pGL4.32  reporter  (Promega, 
Madison,  WI)  which  contains  the  luc2P  reporter  and  hygromycin  resistance.  The  pGL4.32 
vector  was  digested  with  Nhe\  and  Hindlll  (New  England  Biolabs,  Ipswich,  MA)  and  three 
consensus  EREs  spaced  by  three  nucleotides  were  cloned  upstream  of  luc2P  using  the 
following  oligonucleotides:  5'  -CTA  GCG  GTC  ACA  GTG  ACC  TGC  GAG  GTC  ACA 
GTG  ACC  TGC  GAG  GTC  ACA  GTG  ACC  TGC  GA  -  3'  and  5'  -  AGC  TTC  GCA  GGT 
CAC  TGT  GAC  CTC  GCA  GGT  CAC  TGT  GAC  CTC  GCA  GGT  CAC  TGT  GAC  CG  - 
3'.  Successful  cloning  was  verified  by  complete  sequencing  and  the  vector  was  designated 
pGL4.3xERE.  Estrogen  responsiveness  was  validated  by  batch  transfecting  FtEK293  cells 
with  2  ng  of  CMX-ERa  or  CMX-ERp,  45  ng  pGL4.3xERE  vector,  and  40  ng  CMX-P- 
galactosidase  per  well  of  a  48  well  plate.  Cells  were  incubated  24  hr  to  allow  protein 
expression  before  the  addition  of  the  indicated  ligands.  After  24  hr  of  ligand  treatment,  cells 
were  lysed,  firefly  luciferase  substrate  (Promega)  was  added,  and  luminescence  was 
measured  on  a  Victor  X5  microplate  reader  (Perkin  Elmer,  Waltham,  Massachusetts)  using 
luminescence  detection  and  a  700  nm  filter.  To  normalize  data  for  transfection  efficiency,  P- 
galactosidase  expression  was  analyzed  using  the  Tropix  P-galactosidase  detection  kit 
(Applied  Biosystems,  Foster  City,  CA).  Luciferase  counts  were  normalized  to  P-gal  counts 
in  each  well. 

After  characterizing  the  pGL4.3xERE  stable  reporter  vector,  Fls578T-ERa  and  Fls578T-ERp 
cells  were  transfected  with  10  pg  of  the  vector  and  selected  in  125  pg/mL  hygromycin  B  for 
4  weeks.  Individual  colonies  were  selected  using  3  mm  cloning  discs,  expanded,  and 
screened  for  estrogen  induced  luciferase  expression.  One  clone  from  each  cell  line  was 
selected  for  further  characterization,  referred  to  here  as  Fts578T-ERaLuc  and  Fts578T- 
ERpLuc. 
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2.3  Quantitative  western  biots  and  iigand  binding  assays 

For  quantitative  western  blots,  cells  were  split  in  phenol  red  free  DMEM/F12  +  5%  SFS  and 
treated  with  50  ng/mL  Dox  or  vehicle  (water)  24  hr  later.  After  48  hr  treatment,  cells  were 
collected  by  trypsinization,  washed  with  Dulbecco’s  phosphate  buffer  saline  (Invitrogen), 
and  lysed  by  suspension  in  lysis  buffer  (50  mM  Tris  pFl  8.0,  400  mM  NaCl,  10%  glycerol, 
0.5%  triton  X-100,  protease  inhibitors,  and  benzonase).  After  centrifugation,  total  protein 
was  quantified  using  BioRad  Protein  Assay  (BioRad),  and  40  ug  of  protein  was  resolved 
using  SDS-PAGE  and  8%  polyacrylamide  gels.  Proteins  were  transferred  to  a  nitrocellulose 
membrane  for  1.5  hr  at  0.35  A.  Membranes  were  blocked  with  5%  nonfat  milk  and 
incubated  overnight  with  1:1000  anti-FLAG-M2  antibody  (Sigma)  or  1:5000  anti-P-Actin 
(Sigma)  at  4°C.  Membranes  were  then  incubated  with  IRDye  800CW  goat-anti-mouse  IgG 
secondary  antibody  (Licor  Biosciences,  Lincoln,  NE)  for  1  hr  at  room  temperature  and 
visualized  on  a  Licor  Odyssey  near-infrared  gel  reader  (Licor  Biosciences). 

For  ligand  binding  assays,  Fls578T-ERaLuc  and  Fls578T-ERpLuc  cells  were  cultured  in 
phenol  red  free  DMEM/F12  -I-  10%  6x  charcoal  stripped  FBS  (SFS)  for  3  days  prior  to  the 
assay  to  remove  residual  estrogens  from  the  cells.  At  90%  confluence,  cells  were  collected, 
resuspended  in  phenol  red  free  DMEM/F12  -I-  5%  SFS,  and  plated  at  a  density  of  10^  cells/ 
well  on  a  24  well  plate  in  the  presence  or  absence  of  50  ng/mL  Dox.  After  24  hr,  cells  were 
labeled  in  triplicate  with  20  nM  [^Ft]-E2  (89.2  Ci/mmol  specific  activity,  Perkin  Elmer)  in 
the  presence  or  absence  of  450  pM  DES  cold  competitor  for  2  hr  at  37  °C  and  5%  CO2. 
Labeled  cells  were  washed  3  times  with  cold  PBS  +  0.1%  BSA  and  lysed  with  500  pL  SDS 
lysis  buffer  (0.5%  SDS,  0.05  M  Tris-HCl  pH  8.0,  and  1  mM  DTT).  Total  cell  lysate  (400 
pL)  was  mixed  with  5  mL  liquid  scintillation  cocktail  and  [^H]  bound  radioactivity  was 
liquid  scintillation  counted  for  5  min.  Two  additional  wells  of  each  condition  were  used  to 
count  the  cell  number  and  determine  the  total  protein  using  RC  DC  protein  assay  (BioRad, 
Hercules,  CA). 

2.4  Luciferase  assays 

Hs578T-ERaLuc  and  Hs578T-ERpLuc  cells  were  cultured  in  phenol  red  free  DMEM/F12  -1- 
10%  SFS  for  3  days  prior  to  the  assay  to  remove  residual  estrogens  from  the  cells.  Cells 
were  seeded  in  triplicate  at  a  density  of  lO'^  celFwell  on  white  96  well  tissue  culture  plates 
(Fisher)  in  phenol  red  free  DMEM/F12  -I-  5%  SFS  treated  with  50  ng/mL  Dox.  After  24  hr  of 
Dox  treatment,  media  were  replaced  with  treated  media  containing  vehicle  (0. 15%  DMSO) 
or  a  range  of  serially  diluted  ligands.  All  treatments  were  conducted  in  the  presence  and 
absence  of  100  nM  ICI  182,780.  After  treatment  for  24  hr,  cells  were  washed  with  PBS  and 
lysed  with  35  pL  lysis  buffer  (100  mM  K2HPO4,  0.2%  triton  X-100,  pH  7.8).  Lysate  (30  pL) 
was  mixed  1 : 1  with  luciferase  substrate  (Promega)  and  luminescence  was  measured  on  a 
Victor  X5  microplate  reader  (Perkin  Elmer,  Waltham,  Massachusetts)  using  luminescence 
detection  and  a  700  nm  fdter.  Total  protein  (5  pL)  was  quantified  using  BioRad  Protein 
Assay  (BioRad).  EC50  values  were  calculated  using  GraphPad  Prism  Software  (Version 
5.04,  GraphPad  Software  Inc.,  San  Diego,  CA)  and  a  three  parameter  log  versus  response 
nonlinear  regression.  Two  tailed  t-tests  performed  with  GraphPad  Prism  Software  were  used 
to  determine  statistically  significant  differences  from  control  treatments. 

2.5  Gene  expression  analysis 

For  analysis  of  reporter  induction  by  cosmosiin,  Hs578T-ERaLuc  and  Hs578T-ERpLuc 
cells  were  split  in  phenol  red  free  DMEM/F12  -I-  5%  SFS  and  treated  with  50  ng/mL  Dox  for 
48  hr  followed  by  treatment  with  DMSO  (0. 1%),  1  nM  E2,  or  1  pM  cosmosiin  for  4  or  24 
hr.  Total  RNA  was  extracted  using  RNEasy  Plus  Kit  according  to  manufacturer  protocol 
(Qiagen,  Valencia,  CA).  RNA  (2  pg)  was  reverse  transcribed  using  Superscript  II  RT 
according  to  manufacturer  protocol  (Invitrogen),  and  firefly  luciferase  (FLuc)  expression 
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was  determined  by  reverse -transcription  polymerase  chain  reaction  using  primers  shown  in 
Table  1. 

For  quantitative  real-time  PCR  analysis  of  endogenous  target  gene  expression,  Fts578T-ERa 
and  Fts578T-ERp  cells  were  cultured  in  phenol  red  free  DMEM/F12  -I-  10%  SFS  for  3  days 
prior  to  the  assay  to  remove  residual  estrogens  from  the  cells.  Cells  were  split  in  phenol  red 
free  DMEM/F12  -I-  5%  SFS  and  treated  with  50  ng/mL  Dox  for  48  hr  prior  to  ligand 
treatment.  Cells  were  treated  with  Dox  and  ligands  or  vehicle  (0.1%  DMSO)  for  24  hr,  and 
total  RNA  was  extracted  using  RNEasy  Plus  Kit  according  to  manufacturer  protocol 
(Qiagen).  RNA  (2  pg)  was  reverse  transcribed  as  above,  and  quantitative  PCR  was 
performed  using  TaqMan  Prime  Time  custom  designed  assays  (IDT,  Coralville,  lA), 
FastStart  Universal  Probe  Master  Mix  (Roche  Scientific),  and  a  CFX96  instrument 
(BioRad).  Primer  and  probe  sequences  are  shown  in  Table  1.  Data  were  analyzed  using  the 
AACq  method  calculated  by  the  CFX  Manager  Software  (BioRad).  Two  tailed  t-tests 
performed  with  GraphPad  Prism  Software  were  used  to  determine  statistically  significant 
differences  from  control  treatments  using  data  from  three  biological  replicates. 

2.6  Cell  counting  assays 

FIs578T-ERa  and  FIs578T-ERp  cells  were  cultured  in  phenol  red  free  DMEM/F12  -l-  10% 
SFS  for  3  days  prior  to  the  assay  to  remove  residual  estrogens  from  the  cells.  Cells  were 
seeded  at  a  density  of  15,000  cell/well  in  phenol  red  free  DMEM/F12  -l-  5%  SFS  in  triplicate 
in  6  well  tissue  culture  dishes  in  the  presence  or  absence  of  50  ng/mL  Dox.  After  24  hr,  the 
cells  were  treated  with  DMSO  (0.1%)  or  compound  in  the  presence  or  absence  of  50  ng/mL 
Dox.  Media  were  refreshed  every  48  hr,  and  cells  were  counted  after  trypan  blue  exclusion 
using  an  automated  cell  counter  (BioRad)  according  to  manufacturer  protocol. 


3.  Results 

3.1  Generation  of  Hs578T-ERaLuc  and  Hs578T-ERpLuc  reporter  cell  lines 

In  order  to  generate  stable  reporter  breast  cancer  cell  lines,  we  first  cloned  a  construct 
encoding  a  selection  marker  and  a  luciferase  reporter  linked  to  EREs.  The  pGL4.32  vector 
(Promega)  contains  the  luc2P  gene  and  was  modified  to  contain  3  tandem  consensus  EREs 
upstream  of  the  minimal  promoter  (pGL4.3xERE,  Fig.  lA).  Upon  complete  sequencing,  the 
estrogen  responsiveness  of  the  vector  was  validated  in  ER-negative  FIEK293  cells 
transfected  with  full  length  ERa  (Fig.  IB)  or  ERp  (Fig.  1C).  The  pGL4.3xERE  reporter 
showed  extremely  low  background  with  a  65-fold  induction  in  cells  transfected  with  ERa. 
The  ER  antagonist  ICI  182,780  abolished  estrogen  induced  expression,  reducing  the 
luciferase  signal  to  that  of  vehicle  treated  cells.  Cells  transfected  with  ERp  showed  a  15-fold 
induction  of  luciferase  upon  E2  treatment;  ICI  182,780  inhibited  luciferase  expression  in 
both  vehicle  and  estrogen  treated  cells.  The  minimal  background  luciferase  expression  and 
the  selection  marker  conferred  by  the  pGL4.3xERE  vector  made  the  vector  suitable  for 
creating  stable  reporter  cells  lines  for  the  identification  and  characterization  of  ER  selective 
agonists. 

In  order  to  create  stable  ER  reporter  breast  cancer  cell  lines,  an  ER  negative  breast  cancer 
cell  line  engineered  to  express  either  ERa  or  ERp  was  necessary.  Previously,  Secreto  and 
coworkers  created  such  lines  using  FIs578T  cells  [15],  a  triple  negative  breast  cancer  cell 
line  with  a  basal-like  gene  expression  profile  [16].  FIs578T  cells  lack  expression  of  ERa  and 
ERp  providing  a  clean  background  in  which  to  express  ERa  or  ERp.  Using  the  tetracycline 
inducible  system,  two  cell  lines  were  created  in  which  ERa  or  ERp  are  inducibly  expressed 
(Hs578T-ERa  and  Hs578T-ERp  cells,  respectively)  [15].  Hs578T-ERa  and  Hs578T-ERp 
cells  were  transfected  with  the  pGL4.3xERE  vector,  and  individual  clones  were  isolated 
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after  hygromycin  selection.  Over  20  clones  were  screened  for  estrogen  induced  luciferase 
expression  (data  not  shown).  One  clone  from  each  cell  line  was  selected  for  further 
characterization,  referred  to  here  as  Hs578T-ERaLuc  and  Hs578T-ERpLuc.  Additional  ERa 
and  ERp  reporter  clones  were  used  to  verify  reporter  data  obtained  from  Hs578T-ERaLuc 
and  Hs578T-ERpLuc  cells. 

Hs578T-ERaLuc  and  Hs578T-ERpLuc  cells  were  first  characterized  by  assessing  luciferase 
induction  by  ER  ligands  in  the  presence  or  absence  of  the  full  antagonist  ICI  182,780 
(Figure  2).  Cells  were  treated  with  vehicle,  1  nM  E2,  10  nM  DPN  (a  reported  ERp  selective 
agonist),  or  10  nM  PPT  (a  reported  ERa  selective  agonist).  PPT  selectively  activated 
luciferase  expression  in  Hs578T-ERaLuc,  but  DPN  activated  the  reporter  in  both  Hs578T- 
ERaLuc  and  Fls578T-ERpLuc  cells,  though  to  a  lesser  extent  in  Fts578T-ERaLuc  cells.  Co¬ 
treatment  with  ICI  182,780  blocked  luciferase  induction  in  both  cell  lines  (Fig.  2),  and 
luciferase  was  not  induced  in  the  absence  of  Dox  treatment  (data  not  shown). 

Basal  and  E2-induced  luciferase  signals  were  much  higher  in  FIs578T-ERaLuc  cells  when 
compared  to  Hs578T-ERpLuc  cells,  a  trend  observed  in  all  luciferase  assays.  On  average, 
Hs578T-ERpLuc  cells  expressed  630  luciferase  units  per  mg  protein  and  FIs578T-ERaLuc 
expressed  2900  luciferase  units  per  mg  protein  at  saturating  E2  concentrations  (0.1  nM  or 
greater).  A  range  of  luciferase  signals  was  observed  among  the  clones  screened  (data  not 
shown),  suggesting  the  accessibility  of  the  reporter  in  the  chromatin  may  be  responsible  for 
differences  in  luciferase  expression.  In  order  to  verify  Hs578T-ERaLuc  and  FIs578T- 
ERpLuc  cells  had  similar  ER  expression  levels  at  the  Dox  concentration  used  throughout  the 
study  (50  ng/mL),  quantitative  western  blots  were  used  to  compare  ER  expression  in  the 
parent  cell  lines  and  reporter  cell  lines  (Fig.  3A).  Western  blots  with  FLAG  antibody 
demonstrated  similar  ER  expression  in  FIs578T-ERaLuc  and  Hs578T-ERpLuc  cells  and  also 
confirmed  expression  levels  similar  to  the  parent  cell  lines.  In  addition,  whole  cell  ligand 
binding  assays  were  used  to  quantify  the  active  receptor  in  each  cell  line  (Fig.  3B).  ERa 
positive  MCF7  breast  cancer  cells  expressed  -150,000  receptors/cell  which  was  very  similar 
to  reported  values  [17].  Both  FIs578T-ERaLuc  and  Hs578T-ERpLuc  cells  expressed 
-120,000  receptors/cell  after  50  ng/mL  Dox  treatment.  The  comparable  number  of  ERs  per 
cell  suggests  that  differences  in  ER  expression  do  not  account  for  the  higher  luciferase 
signal  observed  Hs578T-ERaLuc  cells.  Higher  luciferase  expression  in  Hs578T-ERaLuc 
cells  may  be  due  to  the  accessibility  of  the  reporter  in  the  chromatin  or  the  enhanced 
transcriptional  activity  of  ERa,  in  agreement  with  previous  findings  that  the  transcriptional 
activity  of  ERa  is  greater  than  that  of  ERp  on  ERE-containing  reporters  [18].  Finally,  the 
reporter  cell  lines  did  not  have  an  altered  morphological  phenotype  compared  to  the  parent 
cell  lines  (Fig.  3C),  and  no  other  phenotypic  changes  due  to  the  integration  of  the  luciferase 
reporter  were  observed  in  Hs578T-ERaLuc  and  Hs578T-ERpLuc  cells. 

3.2  Ligand  selectivity  of  Hs578T-ERaLuc  and  Hs578T-ERpLuc  reporter  cell  lines 

We  next  assessed  ligand  subtype  selectivity  using  these  isogenic  reporter  cell  lines.  All 
luciferase  data  were  normalized  to  the  luciferase  signal  induced  by  a  saturating 
concentration  of  E2  (0. 1  nM)  and  expressed  as  the  percent  transactivation  relative  to  0.1  nM 
E2.  Dose-response  curves  were  obtained  for  E2,  DPN,  and  PPT  to  characterize  the 
sensitivity  of  the  reporter  cells  to  ER  ligands  (Fig.  4).  Cells  were  treated  with  10-fold 
dilutions  of  ligands  and  approximate  EC50  concentrations  for  each  ligand  were  calculated 
from  3  independent  experiments  (Table  2).  The  ratios  of  EC50  values  obtained  from 
Hs578T-ERaLuc  cells  and  Hs578T-ERpLuc  cells  are  also  presented  in  Table  2  and  provide 
a  measure  of  the  selectivity  of  the  ligands.  Higher  o/p  ratios  indicate  selectivity  for  ERp. 

Both  cell  lines  were  highly  sensitive  to  estrogen  (Fig.  4A).  Hs578T-ERaLuc  cells  showed 
EC50  values  near  1  pM;  four  additional  Hs578T-ERaLuc  clones  showed  similar  sensitivities 
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(data  not  shown).  Hs578T-ERpLuc  cells  also  showed  EC50  values  for  estrogen  in  the  pM 
range,  though  the  average  EC50  was  6.5-fold  higher  than  that  of  Hs578T-ERaLuc  cells. 
Similar  differences  in  estrogen  sensitivities  have  been  observed  in  other  ERE-luciferase 
reporter  cell  lines  expressing  ERa  or  ERp  [12-14],  suggesting  the  difference  in  E2 
sensitivity  between  Hs578T-ERaLuc  and  Hs578T-ERpLuc  cells  is  due  to  differences  in  the 
transactivation  of  ERa  and  ERp. 

Next,  dose  responses  to  two  highly  selective  ERa  and  ERp  agonists,  PPT  and  DPN 
respectively,  were  analyzed  using  Hs578T-ERaLuc  and  Hs578T-ERpLuc  cells.  PPT  showed 
nearly  1000-fold  selectivity  for  ERa  (Fig.  4B).  Surprisingly,  PPT  could  activate  reporter 
expression  in  Hs578T-ERpLuc  cells  at  concentrations  greater  than  100  nM,  although  it 
could  not  induce  luciferase  expression  to  the  same  extent  as  E2.  It  has  been  reported  that 
PPT  was  unable  to  induce  an  estrogen  responsive  reporter  in  HEC-1  cells  transfected  with 
ERp  [10]  or  in  HELN-ERa  cells  [12].  DPN  was  not  as  selective  as  PPT  and  could 
maximally  activate  luciferase  expression  Hs578T-ERaLuc  cells  at  100  nM  (Fig.  4C).  DPN 
fully  activated  ERp  at  10  nM.  Though  DPN  has  been  shown  to  have  a  50  to  70-fold  higher 
binding  affinity  for  ERp  [12,  19],  comparison  of  EC50  values  showed  approximately  30-fold 
selectivity  for  ERp  in  these  reporter  assays. 

Next,  the  subtype  selectivity  of  two  natural  phytoestrogens,  liquiritigenin  and  cosmosiin, 
were  analyzed  using  Hs578T-ERaLuc  and  Hs578T-ERpLuc  cells  (Fig.  5).  Liquiritigenin  is  a 
phytoestrogen  derived  from  Glycyrrhizae  umlensis  and  the  most  active  estrogenic 
component  of  MFlOl,  an  herbal  supplement  with  therapeutic  potential  [11].  In  the  initial 
characterization  of  liquiritigenin,  Mersereau  and  coworkers  found  liquiritigenin  showed 
minimal  activation  of  ERa  at  concentrations  up  to  2.5  pM  in  transcriptional  assays  in  U20S, 
HeLa,  or  WAR5  prostate  cancer  cells  transfected  with  ERa  [1 1].  Binding  assays 
demonstrated  that  liquiritigenin  had  a  20-fold  higher  affinity  for  ERp  and  selectivity  was 
proposed  to  be  due  to  selective  recruitment  of  co-activators  to  ERp,  namely  SRC-2  [11]. 
Comparison  of  EC50  values  showed  liquiritigenin  had  a  3.6-fold  selectivity  for  ERp,  and 
maximal  reporter  induction  was  obtained  by  100  nM  liquiritigenin  in  Hs578T-ERpLuc  cells 
and  1  |iM  in  FIs578T-ERaLuc  (Fig.  5A,  Table  2). 

Cosmosiin,  or  apigenin  7-glucoside,  is  a  fiavone  found  in  chamomile  [20]  that  was  identified 
as  an  ER  agonist  that  selectively  induces  ERa/p  and  ERp/p  dimers  as  measured  by 
bioluminescence  resonance  energy  transfer  (BRET)  assays  (unpublished  data).  It  has  a  3- 
fold  higher  binding  affinity  for  ERp  as  measured  by  competitive  ligand  binding  assays  (IC50 
ERa  15.9  pM,  IC50  ERp  3.3  pM,  unpublished  data).  Interestingly,  cosmosiin  induced 
luciferase  expression  to  a  much  greater  extent  than  E2,  an  effect  described  as  supramaximal 
induction  [21].  Even  at  concentrations  up  to  10  pM,  cosmosiin  did  not  saturate  the  luciferase 
output,  and  EC50  values  could  not  be  reasonably  calculated  (Fig.  5B).  Another  FIs578T- 
ERpLuc  clone  treated  with  cosmosiin  also  showed  supramaximal  induction  (data  not 
shown).  Cosmosiin  did  not  induce  luciferase  expression  in  Dox-treated  cells  co-treated  with 
ICI  182,780  or  cells  not  treated  with  Dox  (data  not  shown),  indicating  the  supramaximal 
induction  was  due  to  ERp  activation.  To  determine  if  the  supramaximal  induction  truly 
represented  enhanced  transcriptional  activation,  the  transcript  levels  of  luciferase  were 
assessed  after  4  and  24  hr  treatments  of  E2  and  cosmosiin  (Fig.  5C).  Cosmosiin  did  not 
induce  luciferase  expression  to  a  greater  extent  than  E2  in  either  FIs578T-ERaLuc  or 
FIs578T-ERpLuc  cells,  indicating  alternative  mechanisms  are  responsible  for  the 
supramaximal  effect. 

3.3  Selective  regulation  of  ERa  and  ERp  target  genes  by  ERp  selective  ligands 

We  next  sought  to  validate  the  subtype  selectivity  of  DPN,  PPT,  liquiritigenin  and  cosmosiin 
by  assessing  regulation  of  endogenous  ER  target  genes.  Estrogen  responsive  target  genes  of 
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ERa  and  ERp  were  previously  identified  in  Hs578T-ERa  and  Hs578T-ERp  cells  [15],  and 
two  ERp  target  genes  and  one  ERa  target  gene  were  selected  for  analysis.  Cells  were  treated 
with  50  ng/mL  Dox  for  48  hr  to  induce  expression  of  the  receptors  and  further  treated  with 
the  corresponding  ligands  for  24  hr.  Complement  component  3  (C3,  N]V1_000064)  was  up- 
regulated  in  Hs578T-ERp  cells  upon  E2  treatment  (Fig.  6A).  DPN  and  liquiritigenin  were 
capable  of  inducing  C3  expression  to  a  comparable  level  as  E2  at  concentrations  that  fully 
activate  ERp  with  minimal  ERa  activation,  as  measured  by  reporter  assays  (Fig.  6A). 
Cosmosiin  induced  C3  expression  at  1  pM,  but  not  to  the  same  extent  as  E2,  demonstrating 
cosmosiin  does  not  fully  activate  the  receptor  at  this  concentration.  PPT  slightly  induced  C3 
expression  compared  to  DMSO  in  Hs578T-ERp  cells,  although  PPT  induced  expression  of 
C3  to  a  much  lesser  degree  compared  to  E2.  Repression  of  the  ERp  target  gene  Jagged  1 
{JAGl,  N]Vt_000214)  occurred  to  a  similar  degree  by  E2,  DPN,  liquiritigenin  and  cosmosiin, 
although  100  nM  liquiritigenin  and  1  pM  cosmosiin  do  not  fully  repress  JAGl  expression 
compared  to  E2,  DPN  or  1  pM  liquiritigenin  (Fig.  6B).  Although  the  ERa  selective  agonist 
PPT  slightly  induced  C3  expression  in  Fls578T-ERp  cells,  it  had  no  effect  on  JAGl 
repression,  demonstrating  incomplete  ERp  activation  by  PPT.  To  further  validate  the 
subtype  selectivity  observed  in  reporter  assays,  expression  of  the  ERa  target  gene  alpha-6 
integrin  (ITGA6,  N]V1_000210)  was  determined  after  treatment  of  Hs578T-ERa  cells  with 
E2,  DPN,  PPT,  liquiritigenin  and  cosmosiin.  As  shown  in  Figure  6C,  ITGA6  was  up- 
regulated  by  E2  and  PPT  treatment,  but  DPN  and  liquiritigenin  did  not  fully  activate  its 
expression  at  concentrations  that  showed  selectivity  in  reporter  assays  (10  nM  and  100  nM, 
respectively).  At  1  pM,  liquiritigenin  and  cosmosiin  were  capable  of  activating  ERa,  and 
ITGA6  expression  was  induced  in  Hs578T-ERa  cells. 

Therefore,  the  subtype  selectivity  of  DPN  and  liquiritigenin  observed  in  reporter  cell  lines 
was  validated  by  subtype  selective  regulation  of  endogenous  target  genes.  Cosmosiin, 
however,  activated  expression  of  an  Hs578T-ERa  endogenous  gene  target  at  concentrations 
that  only  slightly  activated  luciferase  reporter  expression  in  Hs578T-ERaLuc  cells. 

3.4  Growth  inhibition  of  Hs578T-ERp  cells  by  liquiritigenin  and  cosmosiin 

We  next  characterized  the  growth  effects  of  liquiritigenin  and  cosmosiin  in  Fts578T-ERa 
and  Fts578T-ERp  cells.  It  was  previously  shown  that  E2  inhibits  the  growth  of  Fts578T-ERp 
cells  [15],  supporting  the  notion  that  the  anti-proliferative  action  of  ERp  may  be  activated  by 
estrogenic  ligands.  We  tested  whether  100  nM  liquiritigenin,  a  concentration  at  which  ERp 
was  selectively  activated,  and  1  pM  cosmosiin  could  also  inhibit  the  growth  of  Fts578T-ERp 
cells.  Fls578T-ERa  and  Fls578T-ERp  cells  were  treated  with  vehicle  (DMSO),  1  nM  E2,  100 
nM  liquiritigenin  or  1  pM  cosmosiin  in  the  presence  or  absence  of  50  ng/mL  Dox  (with  or 
without  ER,  respectively)  for  a  total  of  5  days.  When  ERa  and  ERp  were  not  expressed  (- 
Dox),  the  compounds  had  no  effect  on  the  growth  of  the  cells  (Fig.  7A,  B).  In  contrast,  E2, 
liquiritigenin,  and  cosmosiin  inhibited  the  growth  of  Fts578T-ERp  cells  when  ERp  was 
expressed  (-1-  Dox,  Fig.  7D),  and  there  was  an  approximately  50%  reduction  in  the  number  of 
cells  after  5  days  of  treatment  with  all  three  compounds  (Fig.  7F).  Fts578T-ERa  cells 
showed  slight  inhibition  with  E2  and  liquiritigenin  treatment  when  ERa  was  expressed  (Fig. 
1C),  but  there  was  not  a  statistically  significant  effect  after  5  days  of  treatment  as  measured 
by  2  independent  experiments  (Fig.  7E).  However,  ERa  expression  in  ER  negative  cells 
often  leads  to  growth  inhibition  [22,  23],  and  it  is  likely  that  activation  of  ERa  inhibits  the 
growth  of  Hs578T-ERa  cells.  This  suggests  that  100  nM  liquiritigenin  partially  activates 
ERa  despite  minimal  regulation  of  ITGA6  at  this  concentration. 

4.  Discussion 

ERa  is  an  established  therapeutic  target  for  breast  cancer  treatment,  but  the  development  of 
subtype  selective  estrogenic  ligands  has  gained  interest  with  the  identification  of  ERp  [1]. 
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ERp  opposes  the  actions  of  ERa  suggesting  that  it  may  be  a  potential  therapeutic  target. 
Exogenous  ERp  expression  in  ERa  positive  breast  cancer  cells  impaired  E2  stimulated 
proliferation  [24]  and  tumor  growth  in  xenografts  [25].  In  support  of  the  anti -proliferative 
role  of  ERp,  MCF7  cells  were  more  proliferative  when  ERp  was  knocked  down  [6]. 
Activation  of  ERp  by  subtype  selective  ligands  may  enhance  ERp  growth  repression  without 
stimulating  proliferation  through  ERa;  indeed  ERp  selective  ligands  inhibited  growth  of 
HCl  1  mouse  mammary  cells  [5].  Here,  we  have  also  shown  that  ERp  ligands  can  inhibit  the 
growth  of  breast  cancer  cells  when  ERp  is  expressed.  In  breast  cancer,  however,  ERp 
expression  is  thought  to  decline  during  progression  [26-28]  so  ligands  aimed  at  targeting 
ERp  must  be  highly  selective  and  used  only  in  patients  that  lack  ERa  or  those  with  low 
ERa:ERp  ratios  of  expression.  The  rate  of  ERp  positivity  in  breast  cancer  has  been  reported 
to  range  from  13%  to  83%  [29-32].  In  order  to  effectively  target  ERp  for  cancer  treatment, 
there  is  an  imminent  need  to:  a)  identify  ERp  selective  ligands  with  minimal  side  effects  and 
better  in  vivo  efficacy  and  selectivity,  and  b)  design  clinical  trials  to  recruit  patients  with  low 
ERa:ERp  ratios  in  earlier  stages  of  disease  progression. 

Although  ERp  selective  ligands  have  not  yet  been  used  for  cancer  treatment,  the  therapeutic 
value  of  ERp  has  been  assessed  in  other  diseases.  Two  of  the  most  promising  ERp  selective 
therapies  are  the  ERp  selective  ligand  ERB-041  and  the  herbal  extract  MF-101  [33].  Clinical 
trials  have  been  completed  to  determine  the  efficacy  of  ERB-041  for  treatment  of  Crohn’s 
disease,  endometriosis,  interstitial  cystitis,  and  rheumatoid  arthritis.  Although  results  have 
not  been  published  for  most  of  the  clinical  trials,  results  of  the  rheumatoid  arthritis  trial 
showed  ERB-041  was  well  tolerated  but  did  not  improve  arthritis  symptoms  [34].  MF-101 
also  showed  a  relatively  safe  profile  and  reduced  the  frequency  of  hot  flashes  in  a  phase  II 
clinical  trial  for  treatment  of  post-menopausal  symptoms  [35].  Liquiritigenin  is  the  most 
active  estrogenic  component  of  MF-101[1 1],  suggesting  ERp  selective  ligands  may  prove 
useful  for  treating  post-menopausal  symptoms. 

Strategies  to  identify  ER  subtype  selective  ligands  include  competitive  ligand  binding, 
dimerization,  transcriptional  reporter,  and  proliferation  assays  [21,  36].  Competitive  ligand 
binding  assays  provide  insight  into  binding  affinities  and  are  useful  for  high  throughput 
small  molecule  screening  [37],  but  they  are  limited  because  ligands  can  act  as  agonists  or 
antagonists  and  binding  affinity  does  not  often  reflect  transcriptional  potency.  BRET  assays 
to  measure  receptor  dimerization  have  been  used  to  identify  subtype  selective  ligands  [38], 
but  also  cannot  differentiate  between  agonists  or  antagonists  [39].  Agonists  can  be 
characterized  using  proliferation  assays  in  MCF7  cells,  which  are  highly  sensitive  and 
provide  a  biologically  relevant  endpoint  in  the  context  of  estrogen-sensitive  cells  [40]. 
However,  this  assay  is  limited  by  a  lack  of  specificity,  as  non-estrogenic  mitogens  can 
stimulate  proliferation,  and  cannot  be  used  to  detect  subtype  selective  agonists. 

Transcriptional  assays  can  differentiate  between  agonists  and  antagonists,  overcoming 
limitations  of  binding  and  dimerization  assays.  Mammalian  reporter  cell  lines  useful  for 
identifying  subtype  selective  ligands  have  been  created  from  HeLa  cervical  carcinoma  cells 
[12]  and  HEK293  kidney  cells  [13].  HELN-ERa  and  HELN-ERp  were  generated  from  HeLa 
cells  in  two  steps:  1)  stable  integration  of  ERE-luciferase  to  generate  HELN  cells,  2)  stable 
expression  of  ERa  or  ERp  to  generate  HELN-ERa  and  HELN-ERp  [12].  293/hERa  and  293/ 
hERp  cells  were  generated  by  a  similar  strategy.  Only  one  breast  cancer  reporter  cell  line, 
T47D-KBLuc,  is  available  to  characterize  agonists  in  the  context  of  breast  cancer  cells  [14], 
but  both  ERa  and  ERp  are  expressed,  preventing  identification  of  subtype  selective  ligands. 

In  this  report,  we  described  the  development  of  a  new  set  of  breast  cancer  reporter  cell  lines 
to  characterize  subtype  selective  estrogenic  ligands.  Hs578T-ERaLuc  and  Hs578T-ERpLuc 
cells  were  highly  sensitive  to  E2  with  EC50  values  of  1  pM  and  6.5  pM,  respectively  (Fig. 
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4A).  Similar  E2  sensitivity  was  observed  in  T47D-KBLuc  cells,  which  showed  an 
approximate  EC50  of  3  pM  [14].  Hs578T-ERaLuc  and  Hs578T-ERpLuc  cells  were  more 
sensitive  to  E2  than  HELN-ER  and  293/ER  reporter  cells,  but  all  reporter  cell  lines  showed 
greater  E2  sensitivity  in  ERa  expressing  cells.  HELN-ERa  cells  were  approximately  3  times 
more  sensitive  to  E2  than  HELN-ERp  cells  (EC50  of  0.017  tiM  and  0.068  nM,  respectively) 
[12]  and  293/hERa  cells  were  approximately  4  times  more  sensitive  to  E2  than  cells 
expressing  ERp  (EC50  of  50  pM  and  200  pM,  respectively)  [13].  Although  Hs578T-ERaLuc 
and  Hs578T-ERpLuc  cells  were  not  created  using  the  same  strategy  as  HELN-ER  or  293/ 
hER  reporter  cells  and  likely  have  unique  genomic  integration  of  the  reporter,  similar 
sensitivities  observed  in  all  reporter  cell  lines  suggest  that  this  does  not  inhibit  comparison 
of  subtype  selectivity. 

Reporter  assays  with  two  ER  subtype  selective  ligands  confirmed  that  Hs578T-ERaLuc  and 
Hs578T-ERpLuc  cells  could  be  used  to  differentiate  between  ERa  and  ERp  selective 
ligands.  The  ERp  selective  agonist  DPN  maintained  33-fold  selectivity  in  Hs578T-ERLuc 
cells  (EC50  of  0.26  tiM  for  ERp  and  8.5  nM  for  ERa,  Table  2).  Dose  response  assays  with 
the  ERa  selective  agonist  PPT  revealed  the  sensitivity  of  Hs578T-ERpLuc  cells  (Fig.  4B). 
Although  PPT  was  unable  to  activate  reporter  expression  in  HEC-1  cells  transfected  with 
ERp  [12],  PPT  did  activate  reporter  expression  in  Hs578T-ERpLuc  cells  at  high 
concentrations,  although  not  to  the  full  extent  induced  by  E2.  PPT  reporter  activation  was 
blocked  by  ICI  182,780  co-treatment  (Fig.  2A)  and  did  not  occur  in  the  absence  of  Dox 
treatment  (data  not  shown),  verifying  reporter  activation  was  mediated  by  ERp.  Despite 
activation  of  ERp  at  high  concentrations,  PPT  could  not  fully  activate  reporter  expression  in 
Hs578T-ERpLuc  cells  and  maintained  1000-fold  selectivity  for  ERa. 

Subtype  selectivity  of  two  natural  phytoestrogens,  cosmosiin  and  liquiritigenin,  was  also 
assessed  in  Hs578T-ERaLuc  and  Hs578T-ERpLuc  cells.  Liquiritigenin  maintained 
selectivity  for  ERp  but  to  a  lesser  extent  than  expected,  as  it  has  been  shown  to  minimally 
activate  ERa  in  other  cell  lines  [1 1].  The  discrepancy  in  the  selectivity  of  liquiritigenin  may 
be  due  to  the  enhanced  sensitivity  of  Hs578T-ERaLuc  cells,  differences  in  cofactor 
expression  in  Hs578T  cells,  or  purity  of  the  compound  (our  studies  utilized  commercially 
available  liquiritigenin  and  Mersereau  and  coworkers  [11]  used  extract  from  G.  umlensis). 
The  selectivity  of  cosmosiin  could  not  be  assessed  using  luciferase  assays  due  to 
supramaximal  induction  (Fig.  5B).  Supramaximal  activation  of  estrogen  responsive  reporters 
have  been  described  in  many  systems  [21].  Here,  we  showed  that  supramaximal  induction 
by  cosmosiin  was  not  due  to  enhanced  transcriptional  activation  of  the  reporter  (Fig.  5C). 
Despite  limitations  of  the  reporter  system,  the  subtype  selectivity  of  cosmosiin  could  be 
characterized  by  assessing  target  gene  regulation  in  Hs578T-ERa  and  Hs578T-ERp  cells. 
While  DPN  and  liquiritigenin  maintained  similar  extents  of  selectivity  as  measured  by 
reporter  assays,  cosmosiin  activated  both  ERa  and  ERp  as  measured  by  endogenous  gene 
regulation  (Fig.  6).  Cosmosiin  and  liquiritigenin  induced  similar  growth  inhibitory  effects  as 
E2  in  Hs578T-ERp  cells,  indicating  the  phytoestrogens  could  elicit  ERp  activation  to  a 
similar  extent  as  E2  (Fig.  7). 

Hs578T-ERaLuc  and  Hs578T-ERpLuc  cells  have  several  advantages  for  identifying  ERp 
selective  agonists  in  comparison  to  available  mammalian  reporter  cell  lines.  First,  the 
Hs578T  reporter  cell  lines  have  inducible  expression  of  ERa  and  ERp,  allowing 
determination  of  off-target  reporter  activation  by  assessing  reporter  expression  in  the 
absence  of  Dox.  Second,  Hs578T-ERaLuc  and  Hs578T-ERpLuc  cells  are  highly  sensitive  to 
estrogenic  ligands.  Third,  endogenous  gene  regulation  can  be  used  to  validate  subtype 
selectivity.  Finally,  growth  inhibition  assays  using  Hs578T-ERp  cells  in  the  presence  and 
absence  of  Dox  can  be  used  to  determine  the  biological  endpoint  of  ERp  activation  and 
validate  specificity  of  ligands  to  ensure  they  do  not  have  off-target  cytotoxic  effects.  High 
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throughput  screening  may  be  possible  using  Hs578T-ERaLuc  and  Hs578T-ERpLuc  cells, 
and  luciferase  assay  optimization  using  Hs578T-ERpLuc  cells  has  shown  a  Z  factor  of  0.5 
(data  not  shown),  an  acceptable  range  for  high  throughput  screening  [41].  Therefore, 
Hs578T-ERaLuc  and  Hs578T-ERpLuc  cells  are  useful  for  the  identification  and 
characterization  of  ER  subtype  selective  ligands  that  may  hold  therapeutic  promise. 
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BRET 

bioluminescence  resonance  energy  transfer 

Cos 

cosmosiin 

Dox 

doxycycline 

DPN 

diarylpropionitrile 

E2 

17P-estradiol 

ER 

estrogen  receptor 

PPT 

propyl  pyrazole  triol 

ERE 

estrogen  response  element 

ICI 

ICI  182,780 

Liq 

liquiritigenin 
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Figure  1. 

The  pGL4.3xERE  reporter  construet  is  estrogen  responsive.  (A)  Three  tandem  EREs  were 
inserted  upstream  of  the  lucIP  gene  in  the  pGL4.32  lueiferase  reporter  construct.  HEK293 
cells  were  batch  transfected  with  the  pGL4.3xERE  reporter  construct,  a  P-galactosidase 
construct,  and  full  length  ERa  (B)  or  ERp  (C).  After  allowing  24  hr  for  protein  expression, 
cells  were  treated  in  triplicate  with  vehicle  (DMSO)  or  1  nM  E2  and  vehicle  or  100  nM  ICI 
182,780  (0.15  %  final  DMSO  concentration)  for  an  additional  24  hr.  Raw  lueiferase  units 
(RLUs)  were  normalized  to  P-galactosidase  to  normalize  for  transfection  efficiency.  Error 
bars  represent  standard  deviations. 
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Figure  2. 

ER  subtype  selective  ligands  selectively  induce  luciferase  in  Hs578T-ERaLuc  and  Hs578T- 
ERpLuc  cells.  Hs578T-ERaLuc  (A)  and  Hs578T-ERpLuc  (B)  cells  were  seeded  in  triplicate 
on  96  well  plates  in  the  presence  of  50  ng/mL  Dox  to  induce  ER  expression.  After  24  hr, 
cells  were  treated  with  vehicle  (DMSO),  1  nM  E2,  10  nM  DPN,  or  10  nM  PPT  in  the 
presence  or  absence  of  100  nM  ICI  182,780  (0.15%  final  DMSO  concentration).  Cells  were 
lysed  24  hr  after  ligand  treatment  and  raw  luciferase  units  were  counted.  Error  bars  represent 
standard  deviations.  *  p  values  <  0.05 
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Figure  3. 

Hs578T-ERaLuc  and  Hs578T-ERpLuc  cells  express  similar  levels  of  ER.  (A)  Quantitative 
western  blot  with  Hs578T-ERa  (ERa),  Hs578T-ERaLuc  (ERaLuc),  Hs578T-ERp  (ERP), 
and  Hs578T-ERpLuc  (ERpLuc)  treated  with  vehicle  (-Dox)  or  50  ng/mL  Dox  (-l-Dox).  ER 
expression  was  detected  using  FLAG  antibody  and  quantified  by  normalizing  to  b-actin 
using  the  Licor  Odyssey  near-infrared  gel  reader.  The  normalized  integrated  intensity  for  the 
FLAG  signal  is  shown  below  the  images.  (B)  Ligand  binding  assays  confirmed  the 
quantitative  western  blots.  Fls578T-ERaLuc  and  Fls578T-ERpLuc  cells  were  seeded  in 
triplicate  and  treated  with  vehicle  or  50  ng/mL  Dox  for  24  hr.  Cells  were  labeled  with  20  nM 
[^Ft]-E2  in  the  presence  or  absence  of  cold  competitor  for  2  hr,  washed,  and  total  cell  lysate 
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was  assessed  for  bound  radioaetivity  as  deseribed  in  Materials  and  Methods.  MCF7  cells 
were  included  for  comparison.  Two  additional  wells  of  each  cell  line  and  condition  were 
used  to  determine  the  cell  number  and  the  numbers  of  receptors  per  cell  were  calculated 
based  on  a  1:1  molar  ratio  of  ligand  to  receptor.  The  average  and  standard  deviation  of  three 
independent  experiments  are  shown.  (C)  The  morphology  of  Hs578T-ERaLuc  and  Hs578T- 
ERpLuc  was  similar  to  that  of  the  parent  Hs578T-ERa  and  Hs578T-ERp  cell  lines. 
Representative  phase-contrast  microscopy  images  of  each  cell  line  (lOOX  magnification). 


Biochem  Pharmacol.  Author  manuscript;  available  in  PMC  2012  December  15. 


NIH-PA  Author  Manuscript  NIH-PA  Author  Manuscript  NIH-PA  Author  Manuscript 


Shanle  et  al. 


Page  18 


Hs578T-ERaLuc 
A  •  Hs578T-ER|iLuc 


B 


150i 


100 


o 

m 

(/> 

c 

TO 


50 


10  ‘^0  '-hO  ’^0  ’Oio  »  10  10 


(PPTl  M 

C 


Figure  4. 

Hs578T-ERaLuc  and  Hs578T-ERpLuc  show  subtype  selective  activation.  Dose  response 
curves  of  E2  (A),  PPT  (B),  and  DPN  (C).  Hs578T-ERaLuc  and  Hs578T-ERpLuc  were 
seeded  in  triplicate  and  treated  with  50  ng/mL  Dox  for  24  hr.  Cells  were  then  treated  with  a 
range  of  ligand  concentrations  (0.15%  final  DMSO  concentration)  for  24  hr.  Each  plate 
contained  DMSO,  0.1  nM  E2,  and  100  nM  ICI  182,780  for  controls.  Luciferase  signal  was 
normalized  to  total  protein  in  each  well  and  expressed  as  a  percent  transactivation  relative  to 
signal  obtained  from  saturating  E2  treatment  (0.1  nM).  Each  dose  response  experiment  was 
conducted  at  least  3  times;  data  shown  are  from  one  representative  experiment.  EC50  values 
are  shown  in  Table  2. 
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Figure  5. 

Liquiritigenin  (Liq)  and  cosmosiin  (Cos)  induce  reporter  expression  in  Hs578T-ERaLuc  and 
Hs578T-ERpLuc.  Dose  response  curves  of  liquiritigenin  (A)  and  cosmosiin  (B).  Hs578T- 
ERaLuc  and  Hs578T-ERpLuc  were  seeded  in  triplicate  and  treated  with  50  ng/mL  Dox  for 
24  hr.  Cells  were  then  treated  with  a  range  of  ligand  concentrations  (0.15%  final  DMSO 
concentration)  for  24  hr.  Each  plate  contained  DMSO,  0.1  nM  E2  and  100  nM  ICI  182,780 
for  controls.  Luciferase  signal  was  normalized  to  total  protein  in  each  well  and  expressed  as 
a  percent  transactivation  relative  to  signal  obtained  from  saturating  E2  treatment  (0.1  nM). 
Each  dose  response  experiment  was  conducted  at  least  3  times;  data  shown  are  from  one 
representative  experiment.  EC50  values  are  shown  in  Table  2.  EC50  values  for  cosmosiin 
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could  not  be  determined  because  of  supramaximal  reporter  induction.  The  supramaximal 
induction  by  cosmosiin  was  not  due  to  supramaximal  transcription  of  the  luciferase  reporter 
(C).  Hs578T-ERaLuc  and  Hs578T-ERpLuc  cells  were  treated  with  50  ng/mL  Dox  for  48  hr 
followed  by  treatment  with  DMSO  (0.1%),  1  nM  E2  or  1  pM  cosmosiin  for  4  or  24  hr. 
Firefly  luciferase  (FLuc)  expression  was  determined  by  RT-PCR.  RPL13A  expression  was 
used  to  ensure  equal  loading. 
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Figure  6. 

ERp  selective  ligands  selectively  regulate  ER  target  genes.  Hs578T-ERa  and  Hs578T-ERp 
cells  were  treated  with  50  ng/mL  Dox  for  48  hr  to  induce  ER  expression  followed  by 
treatment  with  the  corresponding  ligands  for  24  hr.  Total  RNA  was  assayed  for  expression 
of  the  ERp  target  genes  C3  and  JAGl  in  Hs578T-ERp  cells  (A,  B  respectively)  and  the  ERa 
target  gene  ITGA6  in  Hs578T-ERa  (C)  cells  by  quantitative  reverse-transcription 
polymerase  chain  reaction.  Target  gene  expression  was  calculated  using  the  AACq  method 
by  normalizing  to  the  ribosomal  protein  RPL13A.  Data  represent  the  average  and  standard 
deviation  of  three  biological  replicates.  *  p  values  <  0.05  compared  to  DMSO  control,  #  p 
values  <  0.05  compared  to  E2  treatment 
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Figure  7. 

Cosmosiin  (Cos)  and  liquiritigenin  (Liq)  inhibit  the  growth  of  Hs578T-ERp  cells.  Hs578T- 
ERa  (A,  C,  E)  and  Hs578T-ERp  cells  (B,  D,  F)  were  seeded  in  6  well  plates  and  treated  with 
vehicle  (A,  B)  or  50  ng/mL  Dox  (C,  D).  After  24  hr,  the  cells  were  treated  with  vehicle 
(0.1%  DMSO)  or  the  indicated  ligands,  and  treatments  were  refreshed  every  48  hr.  Cells 
were  counted  at  the  times  indicated  using  trypan  blue  exclusion.  Comparisons  of  the  cell 
number  on  day  5  are  represented  in  panels  E  (l-ls578T-ERa)  and  F  (Fls578T-ERP).  Data 
represent  two  independent  experiments.  *  p  values  <  0.05 
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Primer  and  Probe  Sequences 


Table  1 


RPL13A 


C3 


JAGl 


ITGA6 


FLuc 


Primerl  5'-TGTTTGACGGCATCCCAC-3' 

Primer  2  5'-CTGTCACTGCCTGGTACTTC-3' 

Probe  5 '  -  CTT  CAG  ACG  CAC  GAC  CTT  GAG  GG  -  3 ' 

Primer  1  5 '  -  AAC  TAG  ATC  ACA  GAGCTGCG-3' 

Primer  2  5'  -  AAG  TCC  TCA  ACG  TTC  CAC  AG  -  3’ 

Probe  5'  -  CGT  TTC  CCG  AAG  TGA  GTT  CCC  AGA  -  3’ 

Primer  1  5'  -  GGA  CTA  TGA  GGG  CAA  GAA  CTG  -  3’ 

Primer  2  5'  -  AAA  TAT  ACC  GCA  CCC  CTT  CAG  -  3' 

Probe  5 '  -  TCA  CAC  CTG  AAA  GAC  CAC  TGC  CG  -  3 ' 

Primer  1  5 '  -  ACC  CGA  GAA  GGA  AAT  CAA  GAC  -  3 ' 

Primer2  5'-CGCCATCTTTTGTGG  GAT  TC  -  3 ' 

Probe  5 '  -  TGG  GTT  GGA  AGG  GCT  GTT  TGT  CA  -  3 ' 
Primer  1  5 '  -  GGC  TGA  ATA  CAA  ACC  ATC  GG  -  3 ' 

Primer  2  5 '  -  CTT  TCT  TGC  TCA  CGA  ATA  CGA  -  3 ' 
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Table  2 

Average  EC50  values  for  ER  Ligands  (M  x  10“^) 


Hs578T-ERaLuc 

Hs578T-ERpLuc 

a/p 

E2 

0.001  (0.0005) 

0.0065  (0.008) 

0.15 

DPN 

8.5(3) 

0.26  (0.02) 

33 

PPT 

0.016(0.001) 

26(21) 

0.001 

Liquiritigenin 

100  (40) 

28(2) 

3.6 
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Polycyclic  aromatic  hydrocarbons  (PAHs)  are  a  diverse  group 
of  widespread  environmental  pollutants,  some  of  which  have  been 
found  to  be  estrogenic  or  antiestrogenic.  Recent  data  have  shown 
that  hydroxylated  PAH  metabolites  may  be  responsible  for  the 
estrogenic  effects  of  some  PAHs.  The  purpose  of  this  study  was 
to  investigate  the  effects  of  several  PAHs,  as  well  as  their  mono- 
hydroxylated  metabolites,  on  estrogen  receptors  (ERs),  ERa  and 
ERp.  Three  parent  PAHs  and  their  monohydroxylated  metabo¬ 
lites  were  each  evaluated  using  transcriptional  reporter  assays  in 
isogenic  stable  cell  lines  to  measure  receptor  activation,  competi¬ 
tive  binding  assays  to  determine  ligand  binding,  and  biolumines¬ 
cence  resonance  energy  transfer  assays  to  assess  dimerization. 
Einaily,  the  estrogenic  effects  of  the  hydroxylated  metabolites  were 
confirmed  by  quantitative  real-time  PCR  of  estrogen-responsive 
target  genes.  Although  the  parent  PAHs  did  not  induce  ERa  or 
ERp  transcriptional  activity,  ail  of  the  monohydroxylated  PAHs 
(1-OH  naphthanol,  9-OH  phenanthrene,  1-OH  pyrene)  selectively 
induced  ERp  transcriptional  activity  at  the  concentrations  tested, 
while  not  activating  ERa.  Additionally,  the  monohydroxylated 
PAHs  were  able  to  competiveiy  bind  ERP,  induce  ERp  homodi¬ 
mers,  and  regulate  ERp  target  genes.  Although  monohydroxylated 
PAHs  appeared  to  have  weak  agonist  activity  to  ERp,  our  results 
showed  that  they  can  elicit  a  biologically  active  response  from  ERp 
in  human  breast  cancer  cells  and  potentially  interfere  with  ERp 
signaling  pathways. 

Key  Words:  polycyclic  aromatic  hydrocarbons;  estrogen  recep¬ 
tors;  monohydroxylated  metabolites;  dimerization;  transcription; 
ligand  binding. 


Polycyclic  aromatic  hydrocarbons  (PAHs)  have  been  of 
increasing  concern  in  the  human  health  field  due  to  their  wide¬ 
spread  dispersion  in  the  environment  and  the  adverse  health 
effects  associated  with  PAH  exposure  (Baird  et  al.,  2005). 
Formed  through  the  incomplete  combustion  of  organic  com¬ 
pounds,  PAHs  can  be  found  in  charbroiled  foods,  cigarette 
smoke,  contaminated  soil,  vehicle  exhaust,  and  in  the  atmosphere 


from  the  by-products  of  industrial  processes.  PAH  exposure  can 
have  several  adverse  effects,  including  carcinogenesis  and  endo¬ 
crine  disruption. 

Although  PAHs  are  a  diverse  group  of  chemicals,  most  are 
metabolized  by  cytochrome  P450s,  a  superfamily  of  enzymes 
that  mediate  the  oxidation  of  lipophilic  substrates  (Anzenbacher, 
2001;  Bauer  et  al.,  1995;  Kim  et  al.,  1998).  The  diol  epoxide 
PAH  metabolites  are  capable  of  inducing  DNA  damage  (Baird 
et  al.,  2005),  and  many  PAHs  have  been  shown  to  be  carcino¬ 
genic  (Bauer  et  al.,  1995;  Kim  et  al.,  1998).  PAHs  can  also  act 
as  endocrine  disrupting  chemicals  by  interfering  with  normal 
estrogen  signaling.  Upon  monohydroxylation,  PAHs  can  induce 
estrogenic  effects  by  directly  interacting  with  estrogen  receptors 
(ERs)  (Arcaro  et  al.,  1999;  Fertuck  et  al.,  2001a,b).  These  data 
suggest  that  the  estrogenic  effects  of  PAHs  are  primarily  medi¬ 
ated  by  the  monohydroxylated  PAH  metabolites. 

ERs,  members  of  the  nuclear  receptor  superfamily  of 
transcription  factors,  exist  in  two  distinct  isoforms,  a  and 
(3.  Encoded  by  separate  genes  on  different  chromosomes, 
ERa  and  ER(3  have  both  overlapping  and  unique  biological 
functions.  The  DNA-binding  domains  share  96%  homology, 
and  ERs  bind  similar  estrogen  response  elements  (EREs) 
to  regulate  transcription  of  target  genes.  The  ligand-binding 
domains  (LBDs),  containing  the  hormone-dependent  activation 
function  (AF-2)  (Tora  et  al.,  1989),  have  55%  identity  and 
have  similar,  but  not  identical,  ligand-binding  pockets  (Pike 
et  al.,  1999).  Upon  ligand  binding,  the  receptors  dimerize  and 
bind  DNA  to  initiate  transcription  of  target  genes  that  mediate 
distinct  biological  effects.  In  the  presence  of  estrogen,  ERa  is 
a  known  driver  of  cell  proliferation,  especially  in  breast  cancer 
cells,  whereas  ER(3  has  been  shown  to  inhibit  ERa-mediated 
cell  proliferation  (Hartman  et  al.,  2006;  Paruthiyil  et  al.,  2004; 
Treeck  et  al.,  2010). 

Given  the  critical  roles  ERs  play  in  regulating  cell  growth  in 
response  to  estrogens,  there  has  been  significant  effort  put  forth 
to  understand  and  predict  the  impacts  of  xenoestrogens  on  ER 
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FIG.  1.  Chemical  stimctures  of  select  polycyclic  aromatic  compounds  and  monohydroxylated  metabolites  studied. 


singaling.  However,  most  studies  have  been  performed  solely 
in  the  context  of  ERa,  with  a  limited  number  of  PAHs  tested. 
Here  we  utilize  several  in  vitro  assays  to  assess  the  effects  of 
three  PAHs  and  their  monohydroxylated  metabolites,  shown  in 
Figure  1,  on  the  transcriptional  activation,  ligand  binding,  and 
dimerization  of  both  ERa  and  ER(3.  Compounds  were  initially 
screened  for  transcriptional  activation  using  a  previously  char¬ 
acterized  pair  of  isogenic  breast  cancer  cell  lines  with  inducible 
expression  of  either  ERa  or  ER(3  and  a  stably  integrated  estro¬ 
gen-responsive  reporter  (Shanle  et  al.,  2011).  These  cell  lines 
provide  a  sensitive  tool  to  directly  compare  the  transcriptional 
induction  of  ERa  and  ER(3.  Next,  bioluminescence  resonance 
energy  transfer  (BRET)  assays  were  performed  to  evaluate  the 
dimerization  status  of  ERs.  BRET  assays  are  able  to  monitor 
protein-protein  interactions  in  a  live,  cell-based  system  (Powell 
and  Xu,  2008;  Tremblay  et  al.,  1999).  Fluorescence  polariza¬ 
tion  experiments  were  utilized  to  generate  competitive  binding 
curves  and  determine  half  maximal  inhibitory  concentration 
(ICjq)  values.  This  provided  a  simple,  yet  specific  way  to  deter¬ 
mine  whether  the  tested  compound  can  compete  with  estrogen 
for  binding  to  ER.  Finally,  compounds  were  evaluated  for  their 
ability  to  upregulate  ER(3  target  genes  via  quantitative  real-time 
PCR  (qPCR). 

Naphthalene,  phenanthrene,  and  pyrene  were  chosen  as  par¬ 
ent  PAH  compounds  for  study  because  they  have  been  detected 


at  high  levels  in  contaminated  environments  (Arcaro  et  al., 
1999),  and  they  are  considered  by  to  be  Priority  Pollutants 
according  to  the  U.S.  Environmental  Protection  Agency.  The 
hydroxylated  metabolites  were  chosen  due  to  their  detection 
after  metabolism  of  the  parent  compound  (Cho  et  al.,  2006; 
Rossbach  et  al.,  2007).  This  is  the  first  study  to  assess  ER  selec¬ 
tive  activity  of  these  PAHs  and  their  hydroxylated  metabolites 
at  the  levels  of  transcriptional  activity  using  isogenic  reporter 
cell  lines,  ligand  binding,  and  dimerization.  The  data  demon¬ 
strate  that  monohydroxylated  PAHs  differentially  interact  with 
ERa  and  ER(3  and  exhibit  stronger  agonistic  activity  toward 
ER(3  compared  with  ERa,  suggesting  that  ER(3-mediated  bio¬ 
logical  processes  need  to  be  evaluated  to  assess  the  outcomes  of 
PAH  exposure  on  humans. 

MATERIALS  AND  METHODS 

Chemicals.  All  PAH  compounds  were  purchased  from  Sigma-Aldrich 
(St  Louis,  MO).  Doxycycline  (Dox)  was  obtained  from  Clontech  (Mountain 
View,  C  A).  ICI 1 82,780  was  obtained  from  Tocris  Bioscience  (Ellisville,  MO). 

Cell  culture  and  reporter  assays.  Cell  culture  media  were  obtained  from 
Invitrogen  (Carlsbad,  CA).  HEK293T  cells  were  maintained  in  Dulbecco’s 
Modified  Eagle’s  Medium  (DMEM)  supplemented  with  10%  Gibco  Fetal 
Bovine  Serum  (FBS;  Invitrogen)  at  37°C  and  5%  CO2.  Hs578T-ERotLuc  and 
Hs578T-ERpLuc  cells  were  previously  created  by  Shanle  et  al.  (201 1)  and  were 


Downloaded  from  http://toxsci.oxfordjournals.org/  at  University  of  North  Carolina  at  Chapel  Hill  on  March  28,  2014 


DIFFERENTIAL  ACTION  OF  PAHs  WITH  ERa  AND  ERp 


361 


TABLE  1 

Primer  and  Probe  Sequences 


RPL13A 

C3 

JAGl 

NRIPl 


Primer  1 
Primer  2 
Probe 
Primer  1 
Primer  2 
Probe 
Primer  1 
Primer  2 
Probe 
Primer  1 
Primer  2 
Probe 


5'-TGT  TTG  ACG  GCA  TCC  CAC-3^ 

5^-CTG  TCA  CTG  CCT  GGT  ACT  TC-3^ 

5'-CTT  CAG  ACG  CAC  GAC  CTT  GAG  GG-3^ 
5'-AAC  TAC  ATC  ACA  GAG  CTG  CG-3^ 

5^-AAG  TCC  TCA  ACG  TTC  CAC  AG-3^ 

5'-CGT  TTC  CCG  AAG  TGA  GTT  CCC  AGA-3' 
5^-GGA  CTA  TGA  GGG  CAA  GAA  CTG-3^ 

5^- AAA  TAT  ACC  GCA  CCC  CTT  CAG-3^ 

5'-TCA  CAC  CTG  AAA  GAC  CAC  TGC  CG-3^ 
5^-AGA  TTC  CCT  GTC  CTC  CTT  CA-3^ 

5^-GGA  AGT  GTT  TGG  ATT  GTG  AGC-3^ 

5'-TGT  GCA  TCT  TCT  GGC  TGT  GTT  TCT  CC-3^ 


maintained  in  DMEM/F12  supplemented  with  L-glutamine  and  10%  Tet-system 
approved  FBS  (Clontech)  at  37°C  and  5%  CO^- 

Reporter  assays  were  performed  as  previously  reported  (Shanle  et  al.,  201 1). 
Briefly,  cells  were  seeded  in  triplicate  at  10'*  cells/well  on  white  96- well  tis¬ 
sue  culture  plates  in  phenol  red-free  DMEM/F12  supplemented  with  5%  char- 
coal-stripped  FBS  treated  with  50ng/ml  Dox.  After  24h,  media  were  removed 
and  replaced  with  media  treated  with  50ng/ml  Dox  and  vehicle  (0.15%  di¬ 
methyl  sulfoxide  [DMSO])  or  PAH  compounds  diluted  in  DMSO.  After  24h 
of  treatment,  the  cells  were  washed  with  30  pi  of  lx  PBS  and  lysed  with  35  pi 
lysis  buffer  (lOOmM  K^HPO^,  0.2%  Triton  X-100,  pH  7.8).  Thirty  microliters 
of  lysate  were  mixed  1:1  with  luciferase  substrate  (Promega,  Madison,  WI), 
and  luminescence  was  measured  with  a  700-nm  filter  on  a  Victor  X5  micro¬ 
plate  reader  (PerkinElmer,  Waltham,  MA).  Total  protein  was  measured  using 
the  Bradford  Method  (Bio-Rad),  and  raw  luciferase  data  were  normalized  to 
total  protein.  Approximate  EC^q  values  were  calculated  using  GraphPad  Prism 
Softwai'e  (Version  5.04;  Graph-Pad  Softwai'e  Inc.,  San  Diego,  CA)  and  a  three- 
parameter  log  versus  response  nonlinear  regression. 

BRET  assays.  The  BRET  assays  were  performed  similarly  to  those  previ¬ 
ously  reported  (Powell  and  Xu,  2008).  Briefly,  HEK293T  cells  were  transfected 
with  BRET  fusion  plasmids  (pCMX-ERa-RLuc  and  pCMX-ERa-YFP  or  pCMX- 
RLuc-ERp  and  pCMX-YFP-ERp).  Twenty-four  hours  after  transfection,  cells 
were  trypsinized  and  resuspended  in  triplicate  in  PBS  at  approximately  50,000  cells 
per  well  in  a  white  96-well  plate.  Cells  were  then  incubated  with  vehicle  (0.6% 
DMSO),  lOnM  E2,  or  monohydroxylated  PAH  compound  for  Ih  at  room  tem¬ 
perature.  Coelenterazine  h  (Promega)  was  added  to  PBS  at  a  final  concentration  of 
5pM.  Emission  measurements  at  460  nm  and  535  nm  were  immediately  taken  on  a 
Victor  X5  microplate  reader  (PerkinElmer).  BRET  ratios  were  calculated  as  previ¬ 
ously  described  (Koterba  and  Rowan,  2006;  Powell  and  Xu,  2008). 

Competitive  binding  assays.  Competitive  binding  assays  were  performed 
using  the  PolarScreen  ERp  Competitive  Binding  Assay  Kit,  Green  (Invitrogen) 
according  to  the  manufacturer’s  protocol.  Recombinant  human  ERp  (20rLM) 
and  fluorescein-labeled  estradiol  were  incubated  for  4  h  with  the  monohydroxy¬ 
lated  PAH  compounds.  Fluorescence  polarization  was  measured  using  a  Victor 
X5  microplate  reader  (PerkinElmer).  Approximate  IC^^  values  were  determined 
by  GraphPad  Prism  Software  (Graph-Pad  Software  Inc.)  from  competitive 
binding  curves. 

Western  blot  analysis.  Western  blots  were  performed  similarly  to  those 
previously  reported  (Shanle  et  al.,  2011)  with  cells  treated  for  48h  with  vehi¬ 
cle  (DMSO)  or  lOpM  monohydroxylated  PAH  compound.  Total  protein  was 
quantified  using  Bio-Rad  Protein  Assay  (Bio-Rad),  35  pg  of  protein  was 
resolved  by  SDS-PAGE,  and  membranes  were  incubated  with  1:1000  anti- 
FLAG-M2  antibody  (Sigma)  overnight  at  4°C.  Membranes  were  then  incu¬ 
bated  with  goat  anti-rabbit  HRP  secondary  antibody  (Licor  Biosciences, 
Lincoln,  NE)  for  1  h  at  room  temperature  and  visualized  using  SuperSignal 
West  Pico  Chemiluminescent  Substrate  (ThermoScientific,  Waltham,  MA) 


on  autoradiography  film.  Membranes  were  then  washed  and  incubated  with 
1:5000  anti-p-Actin  (Sigma)  for  1  h  at  room  temperature,  then  incubated  with 
goat  anti-mouse  HRP  secondary  antibody  (Licor  Biosciences)  for  1  h  at  room 
temperature  and  visualized  using  SuperSignal  West  Pico  Chemiluminescent 
Substrate  (ThermoScientific)  on  autoradiography  film. 

qPCR  analysis.  Hs578T-ERpLuc  cells  were  cultured  in  phenol  red-free 
DMEM/F12  supplemented  with  10%  charcoal-stripped  FBS  for  3  days  prior 
to  experiment  to  remove  any  residual  estrogens.  Cells  were  seeded  into  10-cm 
tissue  culture  plates  in  phenol  red-free  DMEM/F12  supplemented  with  5% 
stripped  serum  and  treated  with  50ng/ml  of  Dox  24  h  prior  to  PAH  treatment. 
Cells  were  then  treated  with  50ng/ml  Dox  plus  0.1%  DMSO  control,  lOnM  E2, 
lOpM  1-OH-naphthalene,  5|iM  9-OH  phenanthrene,  or  5pM  1-OH  pyrene  for 
24  h.  Total  RNA  was  extracted  using  HP  Total  RNA  Kit  (VWR  Scientific,  West 
Chester,  PA)  according  to  the  manufacturer’s  protocol.  One  microgram  of  RNA 
was  reverse  transcribed  using  Superscript  II RT  according  to  the  manufacturer’s 
protocol  (Invitrogen),  and  qPCR  was  performed  using  TaqMan  Prime  Time  cus¬ 
tom  designed  assays  (IDT,  Coralville,  lA),  FastStart  Universal  Probe  Master 
Mix  (Roche  Scientific,  Basel,  Switzerland),  and  a  CFX96  instmment  (Bio-Rad). 
Primer  and  probe  sequences  ai‘e  shown  in  Table  1 . 

Statistical  analyses.  Two-tailed  Student’s  r-tests  were  performed  using 
GraphPad  Prism  version  5.04  for  Windows,  GraphPad  Software  (www. 
graphpad.com). 

RESULTS 

Monohydroxylated  PAHs  Selectively  Activate  ERp  in 

Reporter  Cell  Lines 

In  order  to  test  the  hypothesis  that  hydroxylated  PAHs  may 
have  estrogenic  activity  with  differential  effects  on  ERa  and 
ER(3,  we  first  utilized  Hs578T-ERaLuc  and  Hs578T-ER(3Luc 
reporter  cells  (Shanle  et  al.,  2011).  These  cell  lines  have  induc¬ 
ible  expression  of  ERa  or  ER(3,  respectively,  and  a  stably  inte¬ 
grated  luciferase  reporter  just  downstream  of  three  tandem 
EREs.  Previous  work  has  shown  that  these  cell  lines  are  highly 
sensitive  to  estrogenic  ligands  and  can  be  used  to  distinguish  ER 
subtype  selective  ligands  (Shanle  et  al.,  2011).  In  this  system, 
cells  are  first  treated  with  Dox  to  induce  expression  of  the  recep¬ 
tor,  followed  by  treatment  with  the  corresponding  compounds. 
In  our  initial  experiments  comparing  the  activation  of  ERa  and 
ER(3,  we  observed  that  only  hydroxylated  PAHs  conferred  estro¬ 
genic  activity  at  lOpM  (Fig.  2).  The  monohydroxylated  PAH 
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FIG.  2.  Differential  activation  of  ERa  and  ERp  by  select  monohydroxy- 
lated  PAH  compounds.  (A)  Hs578T-ERaLuc  and  (B)  Hs578T-ERpLuc  stable 
cell  lines  were  treated  in  triplicate  with  lOpM  of  PAH  compound  in  the  pres¬ 
ence  or  absence  of  lOOnM  ICI  182,780  for  24h.  Data  are  expressed  as  fold 
induction  of  raw  luciferase  units  per  mg  protein  over  the  DMSO  control  ±  SD. 
Experiments  were  repeated  at  least  twice.  */?  <  0.01  compared  with  DMSO 
control. 


compounds  were  able  to  induce  the  ERE-luciferase  reporter 
activity  primarily  in  the  Hs578T-ER(3Luc  cells  (Eig.  2B).  In 
these  cells,  1-OH  naphthalene,  9-OH  phenanthrene,  and  1-OH 
pyrene  induced  a  4.2-,  9.7-,  and  8.7-fold  change  over  DMSO 
vehicle  control,  respectively  {p  <  0.01  in  all  cases).  In  contrast, 
only  1-OH  pyrene  induced  the  ERE-luciferase  reporter  activity 
in  the  Hs578T-ERaLuc  cell  line  {p  <  0.01),  but  not  nearly  to 
the  same  degree  as  that  of  17(3-estradiol  (E2)  (Fig.  2A).  The 
ER  antagonist  ICI  182,780  blocked  PAH-induced  expression 
in  all  cases,  reducing  the  luciferase  signal  to  that  of  vehicle- 
treated  cells.  Reporter  expression  induced  by  lOnM  E2  was  not 
fully  blocked  by  ICI  182,780  cotreatment  because  of  the  high 


A. 


Log  Concentration  (M) 

FIG.  3.  Monohydroxylated  PAHs  activate  ERp  in  a  dose-dependent  man¬ 
ner.  Hs578T-ERpLuc  cells  were  treated  with  Dox  for  24h  followed  by  treat¬ 
ment  with  a  range  of  concentrations  of  (A)  1-OH  naphthalene,  (B)  9-OH 
phenanthrene,  or  (C)  1-OH  pyrene.  The  mean  and  SD  shown  are  from  tripli¬ 
cates  of  one  representative  experiment  repeated  twice. 


concentration  and  potency  of  E2.  No  induction  of  reporter  gene 
activity  was  seen  in  control  experiments  in  which  cells  were 
not  treated  with  Dox  (Supplementary  fig.  1),  further  confirming 
ER-mediated  induction  of  the  luciferase  reporter. 
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FIG.  4.  Monohydroxylated  PAH  compounds  selectively  induce  ERp/p  homodimers.  (A)  BRET  data  for  293T  cells  transfected  with  CMX-ERa-RLuc  and 
CMX-ERa-YFP,  showing  no  ERa/a  dimerization  upon  treatment  with  PAH  compounds  in  triplicate.  (B)  BRET  data  for  293T  cells  transfected  with  CMX-RLuc- 
ERp  and  CMX-YFP-ERp,  showing  ERp/p  dimerization  when  treated  with  PAH  compounds  in  triplicate.  The  experiment  was  performed  three  times.  Error  bars 
represent  SEM.  <  0.05  compared  with  DMSO  control. 


We  next  determined  the  dose-dependent  effects  of  the 
hydroxylated  PAHs  in  the  Hs578T-ER(3Luc  cells  (Fig.  3).  The 
half  maximal  effective  concentration  (ECj|,)  values  for  1-OH 
naphthalene  and  1-OH  pyrene  were  found  to  be  approximately 
5.38  and  0.89pM,  respectively.  9-OH  Phenanthrene  proved  to 
be  cytotoxic  at  concentrations  greater  than  lOpM,  and  the 
dose-response  curve  did  not  adequately  saturate;  however,  an 
approximate  ECj|,  value  was  estimated  to  be  >  6.8pM. 

Monohydroxylated  PAHs  Induce  ERfI  Dimers  and  Directly 

Bind  the  Receptor 

To  further  dissect  the  mechanism  through  which  the  monohy¬ 
droxylated  PAHs  activate  ER(3  and  confirm  the  selectivity  of  the 
compounds,  ER  dimerization  induced  by  the  compounds  was 
assessed  using  BRET  assays.  BRET  assays  allow  the  determi¬ 
nation  of  dimer  formation  in  live  cells  by  transfecting  cells  with 
an  energy  donor  (ER  fused  to  Renilla  luciferase)  and  acceptor 
(ER  fused  to  yellow  fluorescent  protein)  (see  Powell  and  Xu, 
2008).  Upon  transfecting  the  cells  with  the  fusion  constructs 
for  ERa  or  ER(3,  9-OH  phenanthrene  and  1-OH  pyrene  were 
shown  to  significantly  induce  ER(3  homodimerization  {p  =  0.02 
and  0.01,  respectively)  (Fig.  4B).  In  contrast,  1-OH  naphthalene 
did  not  significantly  induce  ER(3  dimerization  as  determined  by 
the  BRET  assay  (p  =  0.35).  Following  the  trend  seen  in  the  ERa 
ERE-reporter  assay,  the  monohydroxylated  PAH  compounds 
were  unable  to  induce  ERa  homodimers  (Fig.  4A). 

In  order  to  confirm  that  ER(3  dimerization  and  ERE- 
luciferase  activity  were  directly  induced  by  ligand  binding,  the 
ability  of  the  monohydroxylated  PAH  compounds  to  displace 


fluorescein-labeled  estradiol  from  human  ER(3  was  assessed 
in  a  competitive  binding  assay  (Fig.  5).  The  competition  with 
E2  indicates  that  compounds  directly  bind  to  ER(3  in  the  same 
ligand-binding  pocket  as  E2.  These  competitive  binding  data 
yielded  half  maximal  inhibitory  concentration  (ICj^)  values 
for  9-OH  phenanthrene  and  1-OH  pyrene  at  9.75  and  0.69pM, 
respectively.  In  support  of  the  BRET  results,  1-OH  naphtha¬ 
lene  showed  a  much  lower  affinity  for  ER(3  as  evidenced  by 
Figure  5A,  but  it  was  still  able  to  displace  E2  at  higher  concen¬ 
trations.  The  approximate  ICj^  value  for  1-OH  naphthalene  was 
estimated  at  or  greater  than  0.48pM. 

After  determining  that  monohydroxylated  PAHs  bind  ER(3, 
Western  blots  with  FLAG  antibody  were  used  to  determine  the 
degradation  status  of  the  receptor  (Supplementary  fig.  2),  as  some 
ER  ligands  cause  degradation  of  the  receptor  upon  binding.  These 
Western  blots  confirmed  that  ER(3  was  not  degraded  by  the  mono¬ 
hydroxylated  PAHs  within  48  h  of  treatment. 

Monohydroxylated  PAHs  Exhibit  Estrogenic  Activity  on  ERfl 

Target  Genes 

To  further  validate  the  reporter  assay  and  BRET  assay 
results,  the  regulation  of  endogenous  ER(3  target  genes  was 
assessed.  Estrogen  responsive  target  genes  of  ER(3  were 
previously  identified  in  Hs578T-ER(3  cells  (Secreto  et  al., 
2007).  Two  upregulated  target  genes  (CCS  and  NRIPl)  and 
one  downregulated  target  gene  (JAGl)  were  selected  for 
analysis  by  qPCR  (Fig.  6).  At  lOpM,  1-OH  naphthalene  was 
able  to  induce  CCS  and  NRIPl  expression  2.1-  and  2.2-fold 
over  DMSO,  respectively.  Although  the  increased  expression 
of  CCS  did  not  reach  statistical  significance  (p  =  0.06), 
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NRIPl  was  significantly  upregulated  by  1-OH  naphthalene 
ip  =  0.02).  Treatment  with  5pM  9-OH  phenanthrene  was  able 
to  significantly  induce  CC3  and  NRIPl  expression  2.4-fold 
(p  =  0.02)  and  1.9-fold  (p  <  0.01)  over  DMSO,  respectively. 
Similarly,  5pM  1-OH  pyrene  was  able  to  significantly  induce 
CC3  and  NRIPl  expression  5.6-fold  (p  =  0.02)  and  3.8-fold 
ip  <  0.01)  over  DMSO,  respectively.  Additionally,  all  three 
monohydroxylated  PAH  compounds  were  able  to  downregulate 
the  expression  of  JAGl,  generating  mean  fold  changes  of  0.64 
ip  =  0.02),  0.36  ip  <  0.01),  and  0.32  ip  <  0.01)  over  the  DMSO 
control.  It  is  important  to  note  that  although  all  compounds 
displayed  some  estrogenic  activity  on  the  target  genes  tested, 
the  estrogenic  response  was  not  as  robust  as  that  of  E2. 


DISCUSSION 


FIG.  5.  Monohydroxylated  PAHs  can  bind  ER(3  in  vitro.  Competitive 
binding  curves  for  monohydroxylated  PAH  compounds  displacing  fluorescein- 
labeled  estradiol  from  human  ERjl.  Purified  hERp  and  fluorescein-labeled 
estradiol  were  incubated  for  4h  with  serial  dilutions  in  triplicate  of  (A)  1-OH 
naphthanol,  (B)  9-OH  phenanthrene,  and  (C)  1-OH  pyrene.  Error  bars  repre¬ 
sent  SD. 


Numerous  studies  have  investigated  the  relationship  between 
PAHs,  their  hydroxylated  metabolites,  and  potential  interac¬ 
tions  with  the  ERs,  yet  most  have  focused  on  ERa  (reviewed 
by  Santodonato,  1997).  Hayakawa  etal.  (2007)  reported  estro¬ 
genic  and  antiestrogenic  activity  for  multiple  monohydroxy¬ 
lated  derivatives  of  common  PAHs  in  a  yeast  two-hybrid  assay 
expressing  ERa.  Similar  to  our  findings,  they  also  reported  that 
the  parent  PAH  compounds  lacked  any  estrogenic  or  antiestro¬ 
genic  activity.  Charles  et  al.  (2000)  also  reported  estrogenic 
activity  for  hydroxylated  metabolites  of  the  carcinogen  benzo[a] 
pyrene  (B[a]P)  in  MCP-7  cells,  which  primarily  express  ERa. 
Despite  these  previous  findings,  there  have  been  relatively  few 
studies  comparing  the  effects  of  monohydroxylated  PAHs  on 
the  differential  activation  and  dimerization  of  ERa  and  ER(3. 

Our  results,  consistent  with  prior  studies,  indicate  that  hydroxy¬ 
lated  PAHs  are  the  active  estrogenic  species  and  can  differentially 
activate  either  ERa  or  ER(3.  Although  the  compounds  we  tested 
exhibited  no  interaction  with  ERa,  the  interaction  with  ER(3  is 
novel  and  significant.  Inhibition  of  luciferase  signal  by  the  ER 
antagonist  ICI  182,780,  as  well  as  the  lack  of  luciferase  signal  in 
the  absence  of  Dox,  demonstrates  that  the  results  of  the  reporter 
assay  are  ER(3  mediated.  Competition  with  fluorescein-labeled 
estradiol  indicates  that  these  monohydroxylated  PAH  compounds 
directly  bind  to  ER(3  at  the  same  ligand-binding  pocket  as  E2. 
Eertuck  et  al.  (2001a)  investigated  different  parental  PAH  and 
metabolite  compounds,  and  they  similarly  reported  that  hydroxy¬ 
lated  PAHs  were  able  to  compete  with  estrogen  and  bind  ERs 
with  a  slight  preference  for  ER(3.  Their  data,  consistent  with  our 
findings,  suggest  that  hydroxylated  PAHs  may  preferentially 
affect  ER(3  signaling.  Given  ER(3’s  role  in  normal  development 
and  function  in  reproductive  tissues  as  well  as  in  the  lungs,  colon, 
prostate,  and  cardiovascular  system,  disruption  of  and  interfer¬ 
ence  with  ER(3  signaling  could  have  implications  in  normal  devel¬ 
opment,  as  well  as  in  cancers  and  malfunctions  of  these  tissues. 

In  addition  to  the  reporter  assay  and  competitive  binding  data, 
the  BRET  and  qPCR  data  confirm  that  9-OH  phenanthrene  and 
1  -OH  pyrene  induce  a  biologically  active  ER(3  response  in  this 
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system.  Given  our  data,  1-OH  naphthalene  may  not  necessarily 
induce  ER(3  homodimers  even  at  the  high  concentration  tested 
(lOpM).  In  support  of  these  data,  ligand-hinding  assays  with 
1-OH  naphthalene  demonstrate  a  relatively  low  binding  affinity 
for  ER(3.  Despite  the  negative  BRET  data,  qPCR  for  endogenous 
ER(3  target  genes  suggest  that  1-OH  naphthalene  is  capable  of 
inducing  a  slight  biologically  active  ER(3  response  for  some  ER(3 
target  genes  although  not  to  the  same  extent  as  E2.  Collectively, 
the  data  obtained  for  1  -OH  naphthalene  demonstrate  an  important 
consideration  of  the  in  vitro  assays  used  in  this  study:  different 
assays  have  different  sensitivities  for  detecting  estrogenic 
activity  and  ER  subtype  selectivity.  The  ER(3  homodimerization 
BRET  assay  typically  shows  a  1.5-  to  2-fold  induction  with 
E2  treatment  because  of  high  ligand-independent  dimerization 
(Powell  and  Xu,  2008).  In  addition,  the  BRET  ratios  ultimately 
depend  on  the  conformational  changes  within  the  receptor 
fusion  proteins,  which  allow  for  efficient  energy  transfer,  and 
different  ligands  will  induce  different  conformational  changes, 
thereby  affecting  the  BRET  ratio  output.  Despite  the  lower  fold 
changes  for  the  ER(3  homodimerization  assay,  BRET  has  been 
successfully  used  in  a  high-throughput  manner  to  identify  ER 
dimer  selective  ligands  (Powell  et  ai,  2010)  and,  in  this  study, 
demonstrated  a  signihcant  induction  of  ER(3  homodimerization 
by  two  other  monohydroxylated  PAHs,  1-OH  pyrene  and  9-OH 
phenanthrene. 

Although  each  monohydroxylated  PAH  tested  gave  a  simi¬ 
lar  pattern  of  results,  the  relative  activity  of  each  compound  is 
quite  different.  Our  data  indicate  that  1-OH  naphthalene  is  the 
weakest  ER(3  agonist  among  the  tested  metabolites,  as  demon¬ 
strated  by  low  reporter  gene  output,  a  lack  of  saturation  in  the 
dose-response  reporter  assays,  low  induction  of  ER(3  dimeriza¬ 
tion,  and  a  lower  binding  affinity  for  ER(3.  In  contrast,  1-OH 
pyrene  and  9-OH  phenanthrene  appear  to  be  fairly  efficient 
ER(3  agonists.  Both  ligands  induced  ERE-reporter  gene  activity 
similar  to  E2  and  effectively  displaced  E2  from  the  ER(3  ligand¬ 
binding  pocket.  Both  compounds  also  signihcantly  elicited  ER(3 
homodimerization.  9-OH  Phenanthrene  generated  data  similar 
to  1-OH  pyrene  with  the  exception  that  it  proved  to  be  cytotoxic 
at  concentrations  greater  than  lOpM,  resulting  in  difficulty  to 
obtain  accurate  ECj^,  values.  Despite  the  cytotoxicity  of  9-OH 
phenanthrene  at  high  concentrations,  treatment  with  lower 
concentrations  of  9-OH  phenanthrene  (5|iM)  stimulated  the 
regulation  of  endogenous  ER(3  target  genes  in  Hs578T-ER(3Luc 
cells.  These  data  suggest  that  some  monohydroxylated  PAHs 
can  affect  ER(3-mediated  signaling  prior  to  inducing  general 
cytotoxicity. 


FIG.  6.  Monohydroxylated  PAHs  can  regulate  ERp  target  genes  similar  to 
estradiol.  Expression  of  ERp  target  genes  (CCS,  NRIPl,  and  JAGl)  was  deter¬ 
mined  by  measuring  relative  mRNA  levels  using  qPCR.  RNA  was  collected 
following  treatment  with  0.1%  DMSO,  lOnM  E2,  lOpM  1-OH  naphthalene, 
5pM  9-OH  phenanthrene,  or  5pM  1-OH  pyrene  for  24-h  and  48-h  treatment 
with  50ng/ml  Dox.  Data  are  expressed  as  fold  induction  compared  with  DMSO 
control.  EiTor  bars  represent  SEM.  *p  <  0.05  compared  with  DMSO  control. 


Downloaded  from  http://toxsci.oxfordjournals.org/  at  University  of  North  Carolina  at  Chapel  Hill  on  March  28,  2014 


366 


SEVERS  ETAL. 


Although  our  data  did  not  indicate  that  any  of  the  mono¬ 
hydroxy  lated  PAHs  tested  had  an  effect  on  ERa,  others 
have  reported  ERa  estrogenic  effects  for  these  compounds. 
Hayakawa  et  al.  (2007)  reported  that  all  three  monohydroxy- 
lated  PAHs  exhibited  little  to  no  ERa  estrogenic  activity,  but 
that  1  -OH  pyrene  was  able  to  compete  with  E2  for  ERa  bind¬ 
ing.  Additionally,  Wiele  et  al.  (2004)  reported  that  1-OH  pyrene 
showed  ERa  estrogenic  activity  in  colon  extracts  from  a  simula¬ 
tor  of  the  human  intestinal  microbial  ecosystem.  Discrepancies 
across  these  studies  may  be  due  to  the  use  of  different  assays 
and  cell  lines  to  assess  the  estrogenic  activity. 

Overall,  these  data  suggest  that  common  monohydroxylated 
PAHs  can  interact,  positively  or  negatively,  with  ER  signaling. 
We  can  conclude  from  our  results  and  from  other  studies  that 
hydroxylated  PAHs  are  the  active  estrogenic  species  and  can 
differentially  bind  ERa  or  ER(3,  likely  in  a  cell-  and  tissue- 
specific  manner.  Eew  studies  assessing  the  physiological 
serum  concentrations  of  monohydroxylated  PAHs  have  been 
published,  although  monohydroxylated  PAHs  may  be  used  as 
urine  biomarkers  to  assess  exposure  to  PAHs  (Elovaara  et  al., 
2006).  It  is  therefore  difficult  to  predict  the  concentrations 
of  monohydroxylated  PAHs  that  reach  tissues  such  as  the 
mammary  gland,  and  the  concentrations  shown  to  be  estrogenic 
in  these  studies  may  or  may  not  be  reached  in  the  serum.  Some 
estrogenic  compounds  in  the  diet,  such  as  genistein  found  in  soy 
products,  can  reach  serum  concentrations  near  the  micromolar 
range  (Cassidy  et  al.,  2006).  Ultimately,  the  physiological 
exposure  to  monohydroxylated  PAHs  will  be  a  function  of  both 
exposure  and  metabolic  activity,  which  will  greatly  vary  among 
individuals.  These  in  vitro  studies,  however,  demonstrate  the 
potential  for  monohydroxylated  PAHs  to  impact  ER(3-mediated 
signaling  and  provide  a  framework  for  assessing  the  impacts 
of  other  environmental  chemicals  on  the  dimerization  and 
transcriptional  activities  of  ERa  and  ER(3. 

SUPPLEMENTARY  DATA 

Supplementary  data  are  available  online  at  http://toxsci. 
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Breast  cancers  that  are  negative  for  estrogen  receptor  a  (ERa),  progesterone  receptor,  and  human 
epidermal  growth  factor  receptor  2  are  known  as  triple-negative  breast  cancers  (TNBC).  TNBCs  are 
associated  with  an  overall  poor  prognosis  because  they  lack  expression  of  therapeutic  targets  like 
ERa  and  are  biologically  more  aggressive.  A  second  estrogen  receptor,  ERjS,  has  been  found  to  be 
expressed  in  50%  to  90%  of  ERa-negative  breast  cancers,  and  ERjS  expression  in  TNBCs  has  been 
shown  to  correlate  with  improved  disease-free  survival  and  good  prognosis.  To  elucidate  the  role 
of  ERjS  in  regulating  gene  expression  and  cell  proliferation  in  TNBC  cells,  the  TNBC  cell  line 
MDA-MB-468  was  engineered  with  inducible  expression  of  full-length  ERjS.  In  culture,  ERjS  expres¬ 
sion  inhibited  cell  growth  by  inducing  a  G1  cell  cycle  arrest,  which  was  further  enhanced  by 
17jS-estradiol  treatment.  In  xenografts,  ERjS  expression  also  Inhibited  tumor  formation  and 
growth,  and  17jS-estradiol  treatment  resulted  in  rapid  tumor  regression.  Furthermore,  genomic 
RNA  sequencing  identified  both  ligand-dependent  and  -independent  ERjS  target  genes,  some  of 
which  were  also  regulated  by  ERjS  in  other  TNBC  cell  lines  and  correlated  with  ERjS  expression  in 
a  cohort  of  TNBCs  from  the  Cancer  Genome  Atlas  Network.  ERjS  target  genes  were  enriched  in 
genes  that  regulate  cell  death  and  survival,  cell  movement,  cell  development,  and  growth  and 
proliferation,  as  well  as  genes  involved  in  the  Wnt/jS-catenin  and  the  G1/S  cell  cycle  phase  check¬ 
point  pathways.  In  addition  to  confirming  the  anti-proliferative  effects  of  ERjS  in  TNBC  cells,  these 
data  provide  a  comprehensive  resource  of  ERjS  target  genes  and  suggest  that  ERjS  may  be  targeted 
with  ligands  that  can  stimulate  its  growth  inhibitory  effects.  (Molecular  Endocrinology  27: 
1762-1775,  2013) 


Estrogen  signaling  is  primarily  mediated  by  two  estro¬ 
gen  receptors  (ERs):  ERa  and  ERjS.  ERa  is  expressed 
in  approximately  70%  of  breast  cancers,  and  many  of 
these  cancers  respond  to  endocrine  therapies  that  block 
the  proliferative  action  of  ERa.  Elowever,  approximately 
15%  to  20%  of  all  breast  cancers  lack  expression  of  ERa, 
its  target  gene  progesterone  receptor  (PR),  and  human 
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epidermal  growth  factor  receptor  2  (1-IER2)  and  are  clin¬ 
ically  defined  as  triple-negative  breast  cancers  (TNBCs). 
Eull-length  ERjS  protein  has  been  detected  in  50%  to  90% 
of  ERa-negative  breast  cancers  (1-4),  and  ER/3  expres¬ 
sion  has  been  shown  to  correlate  with  improved  disease- 
free  survival  and  good  prognosis  in  TNBC  (2).  Unlike 
ERa,  for  which  PR  expression  is  indicative  of  receptor 

Abbreviations:  Dox,  doxycycline;  E2,  17/3-estradiol;  ERa,  estrogen  receptor  a;  ER/3,  estro¬ 
gen  receptor  FBS,  fetal  bovine  serum;  GFP,  green  fluorescent  protein;  IPA,  Ingenuity 
pathway  analysis;  HER2,  human  epidermal  growth  factor  receptor  2;  4-OH  Tam,  4-hy- 
droxytamoxifen;  PR,  progesterone  receptor;  RSEM,  RNA  sequencing  expectation  maximi¬ 
zation;  RNA-seq,  RNA  sequencing;  SFS,  charcoal  stripped  FBS;  TCGA,  the  Cancer  Genome 
Atlas;  Tet,  tetracycline;  TNBC,  triple-negative  breast  cancers. 
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expression  and  function,  a  target  gene  or  gene  set  indicat¬ 
ing  ERj8  functionality  has  yet  to  be  identified. 

Like  ERa,  ER/3  is  a  nuclear  receptor  that  regulates 
target  gene  expression  in  estrogen  responsive  tissues,  such 
as  the  mammary  gland.  Multiple  isoforms  of  ERjS  may  be 
expressed  in  the  mammary  gland,  but  the  full-length  re¬ 
ceptor  is  the  only  isoform  able  to  bind  ligand  with  high 
affinity  and  regulate  target  gene  expression  (5,6).  Several 
studies  have  assessed  the  effects  of  full-length  ER/3  expres¬ 
sion  on  the  growth  of  ERa-positive  breast  cancer  cells 
(7-11).  The  results  of  these  studies  demonstrate  that  ERjS 
expression  inhibits  the  proliferative  response  mediated 
by  ERa.  Eewer  reports  have  assessed  the  growth  effects  of 
ER/3  expression  in  breast  cancer  cells  that  lack  ERa,  but 
ER/3  expression  in  ERa-negative  breast  cancer  cells  has 
been  shown  to  inhibit  growth  in  ligand-independent  and 
-dependent  manners  (12-15).  This  leads  to  the  hypothesis 
that  some  ERa-negative  breast  cancers,  including 
TNBCs,  may  benefit  from  therapies  that  target  ERjS  (16). 

In  regard  to  gene  expression,  even  fewer  studies  have 
aimed  to  identify  ERjS  target  genes  in  the  absence  of  ERa. 
Microarray  analyses  comparing  ERa  and  ERjS  target 
genes  in  U20S  osteosarcoma  cells  (17)  and  Hs578T 
breast  cancer  cells  (15)  demonstrated  that  the  two  recep¬ 
tors  have  both  overlapping  and  distinct  target  genes. 
However,  only  ligand-dependent  ER/3  target  genes  have 
been  identified  in  ERa-negative  breast  cancer  cells,  and  a 
comprehensive  assessment  of  both  ligand-independent 
and  -dependent  ERjS  target  genes  in  TNBC  cells  has  not 
yet  been  completed. 

In  an  effort  to  identify  ERjS  target  genes  globally  in 
TNBC  cells  and  assess  the  growth  inhibition  of  ERjS  ex¬ 
pression  in  vitro  and  in  vivo,  we  generated  a  TNBC  cell 
line  with  inducible  expression  of  full-length  ERjS.  We  con¬ 
firmed  ERjS-mediated  growth  inhibition  in  vitro  and  in 
vivo  and  identified  ER/3  target  genes  using  RNA  sequenc¬ 
ing.  We  further  show  that  some  of  the  ERjS  target  genes 
are  associated  with  ERjS  expression  in  a  cohort  of  TNBCs. 
These  data  demonstrate  the  growth  inhibitory  properties 
of  ERjS  and  provide  a  global  view  of  ligand-independent 
and  -dependent  ERjS  target  genes  in  the  absence  of  ERa 
expression.  This  study  also  provides  a  foundation  to  iden¬ 
tify  target  genes  further  that  may  indicate  ERjS  function¬ 
ality  in  TNBCs. 

Materials  and  Methods 

Cell  culture  and  reagents 

Cell  culture  media  were  obtained  from  Invitrogen.  MDA- 
MB-468  cells  were  cultured  in  DMEM  -I-  10%  fetal  bovine 
serum  (FBS;  Invitrogen),  and  BT549  and  HCC1143  cells  were 


cultured  in  RPMI  1640  -I-  10%  FBS.  HCC1143  breast  cancer 
cells  were  obtained  from  the  Leibniz  Institute  DSMZ,  German 
Institute  for  Microorgansims  and  Cell  Culture.  All  cells  were 
cultured  at  37°C  and  5%  CO2.  MDA-MB-468-ERj3  inducible 
cells  were  cultured  in  DMEM  4-10%  Tet-system  approved  FBS 
(Clontech  Mountain  View),  500  mg/L  Zeocin,  and  5  mg/L  Blas- 
ticidin  S  (Research  Products  International).  Before  all  experi¬ 
ments,  MDA-MB-468-ERj8  cells  were  cultured  in  phenol  red 
free  DMEM  -I-  10%  charcoal  stripped  FBS  (SFS)  for  at  least  3 
days  to  remove  residual  estrogens.  Doxycycline  (Dox)  was  ob¬ 
tained  from  Clontech.  The  Flag  and  j8-actin  antibodies  were 
obtained  from  Sigma-Aldrich,  and  the  Hsp90  antibody  was  ob¬ 
tained  from  Santa  Cruz  Biotechnology.  ICI  182,780  and  ERB- 
041  were  purchased  from  Tocris.  All  other  reagents  were  pur¬ 
chased  from  Sigma  Aldrich  unless  otherwise  noted. 

Generation  of  MDA-MB-468-ER/3  inducible  cells 

MDA-MB-468-ERj3  inducible  cells  were  created  as  previ¬ 
ously  described  using  the  TRex  system  (Invitrogen)  (15).  Briefly, 
MDA-MB-468  cells  were  transfected  with  pcDNA6/TR,  and  a 
stable  clone  was  selected  after  selection  with  5  p,g/mL  Blasticidin 
S.  The  MDA-MB-468-Tet  cells  were  then  transfected  with 
pcDNA4/TO-ER/3,  which  encodes  the  530  amino  acids,  full- 
length  ER/3  isoform,  and  selected  with  500  fcg/mL  Zeocin.  Sta¬ 
ble  clones  were  selected  and  screened  for  inducible  expression  of 
ER/3  using  Western  blotting  with  Flag  antibody. 

Ligand  binding  assays  and  Western  blots 

Ligand  binding  assays  were  performed  as  previously  de¬ 
scribed  with  slight  modification  (18).  Briefly,  2  X  10^  MDA- 
MB-468-ER/3  cells  were  seeded  in  a  24-well  plate  and  treated 
with  vehicle  or  50  ng/mL  Dox.  After  48  hours,  cells  were  incu¬ 
bated  with  20  nM  radiolabeled  [^H]-E2  (89.2  Ci/mmol  specific 
activity;  Perkin  Elmer)  in  the  presence  or  absence  of  450  nM 
diethylstilbestrol  cold  competitor.  After  a  2-hour  incubation, 
cells  were  washed  and  lysed  and  the  retained  radioactivity  was 
determined  using  a  liquid  scintillation  counter. 

To  assess  ERjS  expression  using  Western  blots,  MDA-MB- 
468-ER/3  cells  were  treated  with  the  corresponding  Dox  concen¬ 
trations  for  the  given  amount  of  time.  Cells  were  then  washed  in 
PBS  and  lysed  as  previously  described  (18).  After  separation  by 
SDS-PAGE  and  transfer  to  a  nitrocellulose  membrane,  the  mem¬ 
brane  was  incubated  with  the  corresponding  antibody  and  vi¬ 
sualized  by  enhanced  chemiluminescence  on  a  Chemidoc  XRS 
system  (BioRad)  or  by  fluorescence  detection  on  a  Licor  Odys¬ 
sey  gel  reader  (Licor  Biosciences)  as  previously  described  (18). 

Phalloidin  staining 

MDA-MB-468  ERjS#18  cells  were  maintained  in  DMEM 
supplemented  with  10%  stripped  FBS  for  3  days.  The  cells  were 
then  cultured  in  DMEM  4-5%  SFS  with  or  without  Dox  to 
induce  ERjS  expression.  After  24  hours,  the  cells  were  treated 
with  or  without  10  nM  E2.  Forty-eight  hours  before  immuno- 
staining,  the  cells  were  seeded  onto  coverslips  in  a  6-well  plate. 
After  5  days  of  treatment,  the  cells  were  fixed  in  4%  formalde¬ 
hyde  for  15  minutes,  and  the  coverslips  were  washed  three  times 
in  PBS.  Alexa  Fluor  555  Phalloidin  (Cell  Signaling  Technology) 
was  diluted  1:20  (6.6  p,M  stock  concentration  in  methanol)  in 
PBS  and  added  to  the  cells.  After  15  minutes  incubation  at  room 
temperature,  coverslips  were  washed  once  with  PBS.  The  cov- 
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erslips  were  mounted  in  ProLong  Gold  antifade  reagent  with 
4',6-diamidino-2-phenylindole  (Invitrogen).  Fluorescence  was 
detected  using  a  Leica  DM5000B  microscope  (Buffalo  Grove) 
with  the  appropriate  wavelengths. 

Cell  proliferation  and  cell  cycle  assays 

For  cell  proliferation  assays,  1  X  10^  MDA-MB-468-ERj3 
cells  were  seeded  in  triplicate  onto  6-well  plates  in  DMEM  -E 
5%  SES.  After  24  hours,  cells  were  treated  with  vehicle  or  50 
ng/mL  Dox  for  24  hours.  Cells  were  then  treated  with  the 
corresponding  ligands  or  0.1%  DMSO  for  the  given  amount 
of  time.  Media  were  refreshed  every  48  hours.  Cells  were 
trypsinized  and  counted  after  trypan  blue  exclusion  using  an 
automated  cell  counter  (BioRad)  according  to  the  manufac¬ 
turer’s  protocol. 

To  assess  the  cell  cycle  distribution,  3  X  10^  MDA-MB-468- 
ER/3  cells  were  seeded  in  triplicate  onto  6-cm  plates  in  DMEM  -E 
5%  SES.  After  24  hours,  the  cells  were  treated  with  vehicle  or  50 
ng/mL  Dox  for  24  hours  followed  by  treatment  with  0.15% 
DMSO  or  the  corresponding  ligands  for  72  hours.  Cells  were 
collected  by  trypsinization,  fixed  in  cold  95%  ethanol,  and 
washed  in  PBS.  The  fixed  cells  were  then  resuspended  in  pro- 
pidium  iodide  staining  solution  (200  /rg/mL  RNase  A,  50  ju-g/mL 
propidium  iodide,  0.1%  [v/v]  Triton  X-100  in  PBS  -E  1%  BSA) 
and  incubated  overnight  at  4°C.  Samples  were  analyzed  by  flow 
cytometry  at  the  University  of  Wisconsin  Flow  Cytometry 
Laboratory. 

Xenograft  experiments 

All  animal  work  was  performed  in  accordance  with  proto¬ 
cols  approved  by  the  Animal  Care  and  Use  Committee  of  the 
University  of  Wisconsin-Madison.  468-ER/3#18  cells  were  in¬ 
fected  with  retrovirus  generated  from  a  pLNCX-ELuc  (firefly 
luciferase)  vector.  Cells  were  then  cultured  in  regular  media 
containing  400  p-g/mL  G418  (Invitrogen)  for  at  least  eight  pas¬ 
sages  to  generate  luciferase-labeled  468-ERj3#18  cells.  To  assess 
the  effects  of  ER/3  expression  on  tumorigenicity,  2  X  10®  lu¬ 
ciferase-labeled  468-ERj3#18  cells  were  injected  bilaterally  into 
the  inguinal  mammary  fat  pads  of  5-  to  6-week-old  ovariecto- 
mized  athymic  nude-Eoxnlnu  mice  (n  =  4  per  group).  Mice 
were  given  control  diets  or  Dox-containing  diets  (2000  ppm) 
obtained  from  Harlan  Laboratories  immediately  after  the  injec¬ 
tions.  To  assess  tumor  regression  in  response  to  E2,  2  X  10® 
luciferase-labeled  468-ERj3#18  cells  were  injected  as  described 
above,  and  the  tumors  were  allowed  to  form  in  the  absence  of 
Dox  for  10  weeks.  On  day  0,  a  60-day  release  0.05  mg  E2  pellet 
(Innovative  Research  of  America)  was  implanted  into  each 
mouse.  Mice  (n  =  3  per  group)  were  then  separated  randomly 
and  given  a  control  or  Dox-containing  diet.  Luciferase-based 
noninvasive  bioluminescent  imaging  and  analysis  were  per¬ 
formed  as  previously  described  with  an  IVIS  Imaging  System 
(Caliper  Life  Sciences)  (19).  Briefly,  mice  were  anesthetized  and 
injected  ip  with  2  mg  D-luciferin  (10  mg/mL  in  PBS)  (Gold  Bio¬ 
technology).  Imaging  was  completed  between  15  and  20  min¬ 
utes  after  injection.  Eor  bioluminescence  plots,  total  photon  flux 
was  calculated  for  each  mouse  by  using  a  circular  region  of 
interest. 
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Figure  1.  MDA-MB-468-ERj3  cells  express  ERjS  after  Dox  treatment. 
Dox-dependent  expression  of  ER^  was  Initially  characterized  In  two 
MDA-MB-468-ER/3  clones  (468-ER/3#18  and  468-ER/3#32).  (A)  A 
Western  blot  using  Flag  antibody  shows  that  both  clones  express 
detectable  levels  of  ERj3  after  treatment  with  as  little  as  5  ng/mL  Dox 
for  48  hours.  Further  In  vitro  studies  were  performed  with  50  ng/mL 
Dox.  (B)  A  Western  blot  using  Flag  antibody  shows  that  treatment  with 
50  ng/mL  Dox  for  as  little  as  4  hours  induces  expression  of  ERjS.  (C) 
Ligand  binding  assays  performed  with  radiolabeled  E2  demonstrate 
that  both  clones  express  similar  levels  of  ERj8  corresponding  to 
approximately  90  000  receptors/cell  after  24  hours  of  Dox  treatment. 
The  data  represent  the  mean  and  SD. 


RNA  sequencing  (RNA-seq)  of  ER/3  target  genes 
and  pathway  analysis 

To  identify  ERjS  target  genes  globally,  468-ERj3#18  cells 
were  treated  with  vehicle  or  50  ng/mL  Dox  for  48  hours  fol¬ 
lowed  by  treatment  with  0.1%  DMSO  or  10  nM  E2  for  4  hours. 
Total  RNA  was  extracted  using  a  RNEasy  Plus  Kit  according  to 
manufacturer  protocol  (Qiagen),  and  three  independent  exper¬ 
iments  were  performed.  The  total  RNA  was  submitted  to  the 
University  of  Wisconsin  Biotechnology  Center  for  RNA  quality 
analysis,  sequencing  library  generation,  and  sequencing.  An  II- 
lumina  HiSeq  2000  was  used  to  generate  100-bp  single-end 
reads.  Reads  were  trimmed  for  quality  (15  bp  from  the  5'  end 
and  2  bp  from  the  3'  end)  and  aligned  to  the  University  of 
California,  Santa  Cruz  build  hgl9  genome  using  RNA  sequenc¬ 
ing  expectation  maximization  (RSEM).  It  has  now  been  well 
recognized  that  discarding  RNA-seq  reads  that  align  to  multiple 
genes  introduces  bias  in  expression  estimates  of  gene  families 
and  genes  containing  or  inside  of  repetitive  elements  (20-23). 
RSEM  implements  a  principled  and  biologically  motivated 
probabilistic  model  for  this  problem  (20,  24).  It  is  based  on  a 
generative  model  of  the  RNA-seq  protocol  that  takes  into  ac¬ 
count  aspects  such  as  sequencing  error,  fragment  length  distri¬ 
butions,  and  nonuniformities  in  the  distribution  of  read  loca- 
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tions  along  transcripts.  The  RSEM  alignment  results  are 
presented  in  Supplemental  Table  2,  published  on  The  Endocrine 
Society’s  Journals  Online  web  site  at  http://mend.endojournals. 
org.  Eor  the  identification  of  differentially  expressed  genes,  the 
fractional  counts  were  rounded  to  the  nearest  integer,  and 
DESeq  (25)  was  used  to  assess  differentially  expressed  genes  at  a 
1.5-fold  change  cutoff  and  a  false  discovery  rate  <0.01.  The 
filtered  gene  lists  for  each  comparison  are  provided  in  Supple¬ 
mental  Materials  and  Methods.  The  raw  data  will  also  be  sub¬ 
mitted  to  the  Gene  Expression  Omnibus.  Pathway  enrichments 
and  functional  analyses  of  the  gene  lists  were  determined  using 
ingenuity  pathway  analysis  (IPA)  software  (Ingenuity  Systems, 
www.ingenuity.com).  The  DESeq  and  IPA  results  for  all  com¬ 


parisons  are  available  in  the  Supplemental  Materials  and 
Methods. 

Quantitative  PCR 

468-ER/3#18  cells  were  treated  with  DMSO,  10  nM  E2,  or 
100  nM  ERB-041  as  described  for  the  initial  RNA-seq  experi¬ 
ment.  Eor  validation  in  HCC1143  cells,  the  cells  were  infected 
with  retrovirus  generated  from  pLNCX-GEP  or  pLPCX-ERj3 
vectors.  Twenty-four  hours  after  infection,  the  media  were  re¬ 
placed  with  phenol  red  free  RPMI  1640  -I-  5%  SES.  After  an 
additional  24  hours,  cells  were  treated  with  0.1%  DMSO  or  10 
nM  E2  for  4  hours.  Total  RNA  was  extracted  as  described  above 
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Figure  2.  ERj3  expression  and  E2  treatment  inhibit  the  growth  of  MDA-MB-468-ERJ3  cells.  (A)  Two  MDA-MB-468-ERJ3  clones  (#18  and  #32)  were 
treated  in  triplicate  with  vehicle  or  Dox  for  24  hours  followed  by  treatment  with  DMSO  or  1  nM  E2.  The  total  number  of  viable  ceils  was 
determined  using  trypan  blue  exclusion  after  the  given  number  of  days.  * ,  P  <  .01  compared  to  DMSO  control;  #,  P  <  .01  compared  to 
DMSO-EDox  control.  (B)  Treatment  with  ER  antagonists  rescue  the  growth  inhibition  mediated  by  ER^.  468-ERj3#1 8  cells  were  treated  in  triplicate 
with  vehicle  or  Dox  for  24  hours  followed  by  treatment  with  DMSO  or  E2  in  the  presence  or  absence  of  1 00  nM  ICI  or  4-OH  Tam.  The  total 
number  of  viable  cells  was  determined  by  trypan  blue  exclusion  after  5  days.  *,  P  <  .01  compared  to  DMSO  control;  #,  P  <  .01  compared  to 
DMSO-EDox  control.  (C)  ERjS  expression  and  E2  treatment  alter  the  morphology  of  the  468-ERp#1 8  cells.  Cells  were  treated  with  vehicle  or  Dox 
for  24  hours  followed  by  treatment  with  DMSO  or  E2  for  5  days  and  subsequently  stained  with  Alexa  Fluor  B5B  Phalloidin  {red),  which  stains  actin 
filaments  to  highlight  the  cytoplasm,  and  4',6-diamidino-2-phenylindole  {blue)  for  nuclear  staining.  (D)  ER/3  expression  and  E2  treatment  induce  a 
G1  cell  cycle  arrest.  468-ER/3#1 8  cells  were  treated  in  triplicate  with  vehicle  or  Dox  for  24  hours  followed  by  treatment  with  DMSO  or  E2  in  the 
presence  or  absence  of  100  nM  ICI  for  72  hours.  The  proportion  of  cells  in  each  phase  of  the  cell  cycle  was  determined  using  flow  cytometry  as 
described  in  the  Materials  and  Methods.  For  the  proportion  of  cells  in  G1 ;  *,  P  <  .01  compared  to  DMSO  control;  #,  P  <  .01  compared  to 
DMSO-EDox  control.  The  data  represent  the  mean  and  SD. 
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Figure  3.  ERjS  expression  inhibits  tumor  growth  in  xenografts.  (A)  468-ER/S#18  cells  were 
labeled  with  luciferase  and  injected  into  the  mammary  fat  pads  of  nude  mice  (n  =  4  per  group). 
Mice  were  then  provided  a  control  or  Dox-containing  diet.  The  tumor  growth  was  monitored 
using  bioluminescence  imaging  as  described  in  the  Materials  and  Methods.  (B)  Quantitative  plot 
of  bioluminescence  for  each  group  over  time.  (C)  Final  tumor  volumes  for  each  group.  The  data 
represent  the  mean  and  SEM  of  the  tumors  from  each  group.  * ,  P  <  .05  compared  to  control. 


and  reverse-transcribed  using  Superscript  II RT  according  to  the 
manufacturer’s  protocol  (Invitrogen).  Quantitative  PCR  was 
performed  using  TaqMan  Prime  Time  custom  designed  assays 
(IDT),  FastStart  Universal  Probe  Master  Mix  (Roche  Scientific), 
and  a  CFX96  instrument  (BioRad).  Primer  and  probe  sequences 
are  shown  in  Supplemental  Table  1.  Data  were  analyzed  using 
the  AACq  method  calculated  by  CFX  Manager  Software  (Bio- 
Rad).  Student  t  tests  were  performed  with  GraphPad  Prism  Soft¬ 
ware  (Version  5.04;  GraphPad  Software  Inc). 

Correlation  analysis  of  ER/3  target  genes  and  ER/3 
expression  in  human  tumor  samples 

To  compare  the  RNA  sequencing  results  further  with  previ¬ 
ous  studies  of  TNBC,  the  Agilent  mRNA  expression  microarray 
data  from  the  Cancer  Genome  Atlas  (TCGA)  database,  among 
which  n  =  43  samples  were  identified  as  TNBC,  was  examined 
to  assess  any  correlations  between  the  expression  of  ERjS  and  its 
target  genes  (26).  Using  all  of  the  genes  identified  in  the  com¬ 
parison  of  468-ERj3#18  cells  treated  with  E2  versus  E2-FDox 
(comparison  IV,  Eigure  5A)  with  false  discovery  rate  <0.05,  a 
total  of  2472  genes  were  present  in  both  the  microarray  data  and 
our  gene  list.  The  correlation  coefficient  was  calculated  for  the 
expression  of  each  of  the  genes  and  that  of  ERjS  based  on  the 
microarray  data.  To  assess  the  strength  of  the  correlation,  the  P 
values  were  calculated  for  a  two-sided  alternative  that  these 
correlation  coefficients  are  significantly  different  from  zero 
(meaning  no  correlation).  To  account  for  multiple  testing,  these 
P  values  were  then  converted  to  the  so-called  q  values  (27,  28). 
Eor  those  genes  with  q  <  .05  for  the  correlation  and  a  fold 
change  >1.5  in  the  RNA  sequencing  results,  the  correlation 
coefficient  determined  for  the  array  data  was  graphed  against 
the  fold  change  observed  for  the  gene  in  the  468-ERj8#18  cells 
treated  with  E2  versus  E2-FDox. 


Results 

Inducible  expression  of  ER/3  and 
E2  treatment  induces  a  G1  cell 
cycle  arrest  and  inhibits  the 
growth  of  MDA-MB-468  cells  in 
vitro  and  in  vivo 

To  assess  the  growth  and  gene  ex¬ 
pression  effects  of  ERjS  expression 
and  activation  in  TNBC  cells,  MDA- 
MB-468  cells  were  engineered  with 
inducible  expression  of  full-length 
ERjS  using  the  TRex  system  from  In¬ 
vitrogen.  After  creating  MDA-MB- 
468-Tet  cells,  which  express  the  tet¬ 
racycline  (Tet)  repressor,  these  cells 
were  transfected  with  a  construct 
encoding  Tet-operated  full-length 
Elag-tagged  ERjS.  Two  MDA-MB- 
468-ERj8  clones  with  inducible  ERjS 
expression  were  selected  for  further 
characterization  (486-ERj8#18  and 
468-ERj8#32).  Eirst,  ERjS  expres¬ 
sion  was  assessed  in  response  to  increasing  Dox  concen¬ 
trations  (Eigure  lA).  As  little  as  5  ng/mL  Dox  treatment 
for  48  hours  induced  expression  of  the  receptor,  and  sub¬ 
sequent  in  vitro  experiments  were  performed  with  50 
ng/mL  Dox.  Next,  a  time  course  was  established  for  the 
Dox-inducible  expression  of  ERjS  (Eigure  IB).  ERjS  was 
maximally  expressed  after  24  hours,  and  this  expression 
level  was  sustained  after  48  hour  and  72  hour  treatments. 
Einally,  the  expression  level  was  quantified  using  radio- 
labeled  ligand  binding  assays  (Eigure  1C).  After  Dox 
treatment,  both  clones  expressed  approximately  90  000 
receptors  per  cell,  which  is  comparable  to  MCE7  breast 
cancer  cells  that  express  ~150  000  ERa  molecules  per 
cell  (18). 

After  confirming  Dox-inducible  expression  of  ERjS  in 
486-ERj3#18  and  468-ERj3#32  cells,  we  assessed  the 
growth  effects  of  ERjS  expression  and  activation  by  17jS- 
estradiol  (E2)  (Eigure  2A).  The  number  of  viable  cells  was 
significantly  inhibited  by  ERjS  expression  and  E2  treat¬ 
ment  for  both  clones.  468-ERjS#18  showed  slight  growth 
inhibition  when  ERjS  was  expressed  in  the  absence  of 
ligand.  Dox  and  E2  treatment  had  no  effect  on  the  growth 
of  MDA-MB-468-Tet  cells  (Supplemental  Eigure  1).  To 
confirm  that  the  growth  inhibitory  effects  observed  in 
MDA-MB-468-ERjS  cells  were  mediated  by  ERjS,  468- 
ERjS#18  cells  were  treated  with  vehicle  or  Dox  and 
DMSO  or  E2  in  the  presence  or  absence  of  the  antagonist 
ICI  182,780  (ICI)  or  the  selective  estrogen  receptor  mod¬ 
ulator  4-hydroxytamoxifen  (4-OH  Tam)  (Eigure  2B). 


The  Endocrine  Society.  Downloaded  from  press.endocrine.org  by  [SlindividualUser.displayName}]  on  28  March  2014.  at  12:20  For  personal  use  only.  No  other  uses  without  permission. .  All  rights  reserved. 


doi:  10.1210/me.2013-1164 


mend.endojournals.org  1767 


Dox  treatment  alone  significantly  inhibited  the  growth  of 
468-ERj3#18  cells,  and  this  effect  was  rescued  by  ICI  or 
4-OH  Tam  co-treatment.  Treatment  with  E2  further  in¬ 
hibited  the  growth  of  the  cells,  and  this  inhibition  was 
also  rescued  by  ICI  or  4-OEI  Tam  co-treatment.  The  ERjS- 
selective  ligand,  ERB-041  (29),  also  inhibited  cell  growth 
similar  to  E2  (Supplemental  Eigure  2).  Interestingly,  the 
morphology  of  the  cells  was  altered  after  5  days  of  treat¬ 
ment  with  Dox  and  E2  (Eigure  2C),  suggesting  that  the 
cells  may  have  undergone  some  type  of  cell  cycle  arrest. 
Indeed,  an  analysis  of  the  cell  cycle  phase  distribution  by 
flow  cytometry  revealed  an  increase  in  the  proportion  of 
cells  in  the  G1  phase  after  Dox  treatment  alone,  which 
was  further  increased  by  E2  treatment  (Eigure  2D).  ER/3 
expression  and  ligand  treatment  did  not  induce  detectable 
levels  of  apoptosis  (data  not  shown). 

Next,  we  aimed  to  assess  the  effects  of  ER/3  expression 
and  activation  on  tumor  growth  of  MDA-MB-468  cells. 
Eirst,  soft  agar  colony  formation  assays  revealed  that  ER/3 
expression  alone  completely  blocked  the  formation  of 
colonies  (Supplemental  Eigure  3).  Xenograft  experiments 
were  then  performed  using  468-ER/3#18  cells  retrovirally 
labeled  with  luciferase  to  monitor  tumor  growth  by  bio¬ 
luminescence  imaging  (Eigure  3).  After  injecting  cells  into 
the  mammary  fat  pads,  the  mice  were  given  a  control  diet 
or  a  Dox-containing  diet,  and  the  tumor  size  was  moni¬ 
tored  using  bioluminescence  imaging  (Eigure  3A).  In  the 
presence  of  Dox,  the  growth  of  the  tumors  was  reduced 
and  one  mouse  did  not  develop  tumors  (Eigure  3,  A  and 
B).  The  final  tumor  volume  was  significantly  smaller 
when  ER/3  was  expressed  (Eigure  3C),  and  Dox  treatment 
had  no  effect  on  the  growth  of  control  MDA-MB-468-Tet 
tumors  (Supplemental  Eigure  4).  To  assess  the  effects  of 
ER/3  expression  and  activation  on  tumor  growth,  lu- 
ciferase-labeled  468-ER/3#18  cells  were  allowed  to  form 
tumors  for  10  weeks.  Mice  were  then  implanted  with  a 
slow  release  E2  pellet  and  given  a  control  or  Dox-contain- 
ing  diet.  The  tumors  rapidly  regressed  only  when  ER/3 
was  expressed  (Eigure  4).  These  data  demonstrate  the 
growth  inhibitory  effects  of  ERj8  expression  and  activa¬ 
tion  in  TNBC  cells  in  vivo. 

ER/3  target  gene  identification  using  RNA-seq 

The  MDA-MB-468-ERj8  inducible  cells  are  useful  for 
identifying  the  ligand-independent  and  ligand-dependent 
target  genes  of  ER/3  in  the  absence  of  ERa.  To  identify 
direct  ER/3  target  genes,  468-ERj8#18  cells  were  treated 
with  or  without  Dox  followed  by  treatment  with  E2  for  4 
hours.  RNA-seq  was  subsequently  performed  as  de¬ 
scribed  in  the  Materials  and  Methods.  Reads  were  aligned 
using  RSEM  (20),  and  the  results  of  the  alignment  are 
shown  in  Supplemental  Table  2.  Differentially  expressed 
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Figure  4.  Activation  of  ERj3  in  xenografts  causes  tumor  regression. 

(A)  468-ERj3#18  cells  were  labeled  with  luciferase  and  injected  into  the 
mammary  fat  pads  of  nude  mice.  After  allowing  tumors  to  form  for  10 
weeks,  slow  release  E2  pellets  were  implanted  into  each  mouse.  Mice 
were  separated  into  two  groups  (n  =  3  per  group)  and  provided  a 
control  or  Dox-contalning  diet.  The  tumor  growth  was  monitored 
using  blolumlnescence  Imaging  and  caliper  measurements  as  described 
in  the  Materials  and  Methods.  (B)  Quantitative  plot  of  bioluminescence 
for  each  group  over  time.  The  data  represent  the  mean  and  SEM  of 
the  tumors  from  each  group.  (C)  Quantitative  plot  of  tumor  size  as 
determined  by  caliper  measurements  for  each  group  over  time.  The 
data  represent  the  mean  and  SEM  of  the  tumors  from  each  group. 

* ,  P  <  .05  compared  to  control. 

genes  were  determined  for  the  four  conditions  by  DESeq 
(25)  (Eigure  5A,  1.5-fold  change  cutoff).  No  genes  were 
differentially  expressed  after  E2  treatment  in  the  absence 
of  ER/3  (DMSO  vs  E2,  comparison  I),  confirming  that 
MDA-MB-468  cells  are  not  estrogen  responsive.  Only 
109  genes  were  differentially  expressed  after  Dox  treat¬ 
ment  (DMSO  vs  DMSO -E ER/3,  comparison  II),  whereas 
481  genes  were  differentially  expressed  after  E2  treatment 
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Figure  5.  Global  identification  of  ER/3  target  genes  using  RNA-seq.  468-ERj8#18  cells  were  treated  with  vehicle  or  Dox  for  48  hours  followed  by 
treatment  with  DMSO  or  1 0  nM  E2  for  4  hours.  Total  RNA  was  extracted  and  sequenced  as  described  in  the  Materials  and  Methods.  Differentially 
expressed  genes  were  identified  using  DESeq  with  a  1 .5-fold  change  cutoff  and  a  false  discovery  rate  <0.01 .  (A)  A  matrix  describing  the  four 
contrasts  performed  to  identify  differentially  expressed  genes.  No  genes  were  differentially  expressed  in  the  DMSO  vs  E2  group  (I),  109  were 
differentially  expressed  in  the  comparison  between  DMSO  vs  DMSO  +  ERjS  (II),  481  genes  were  differentially  expressed  between  DMSO  +  ERjS  vs 
E2-I-ER^  (III),  and  930  genes  were  differentially  expressed  between  E2  vs  E2-l-ERj3  (IV).  (B)  A  Venn  diagram  showing  the  overlap  among  genes 
identified  in  the  three  groups.  (C)  A  Venn  diagram  showing  the  overlap  between  the  up-regulated  genes  identified  in  DMSO  vs  DMSO  +  ERjS  and 
DMSO-i-ER/3  vs  E2-I-ER/3.  Most  ERj8  target  genes  were  up-regulated,  and  approximately  half  of  the  DMSO  vs  DMSO-l-ER/3  genes  were  identified 
in  DMSO-EERjS  vs  E2-EERJ3.  (D)  A  Venn  diagram  showing  the  overlap  between  the  down-regulated  genes  identified  in  DMSO  vs  DMSO-l-ER/3  and 
DMSO-)-ERj3  vs  E2-l-ERj3.  Eewer  ERjS  target  genes  were  down-regulated,  and  only  five  genes  were  identified  in  both  groups. 


in  the  presence  of  ERjS  (DMSO+ERjS  vs  E2+ERj3,  com¬ 
parison  III),  suggesting  that  most  ERjS  target  genes  are 
ligand-dependent  in  these  cells.  The  top  genes  regulated  in 
comparisons  II  and  III  are  shown  in  Tables  1  and  2,  re¬ 
spectively.  A  comparison  of  E2  treatment  with  or  without 
Dox  (E2  vs  E2+ERj8,  comparison  IV)  revealed  930  dif¬ 
ferentially  expressed  genes,  potentially  capturing  the 
total  change  in  gene  expression  in  response  to  both 
ERjS  expression  and  activation  by  E2.  The  top  genes  in 
this  comparison  are  shown  in  Table  3. 

Most  of  the  genes  in  comparisons  II  and  III  were  cap¬ 
tured  in  comparison  IV  (Eigure  5B).  In  addition,  most 
genes  identified  in  comparisons  II  and  III  were  up-regu- 
lated,  and  about  half  of  the  potentially  ligand-indepen¬ 
dent  genes  (comparison  II)  were  also  regulated  in  a  ligand- 
dependent  manner  (Eigure  5C).  Eewer  genes  were  found 
to  be  down-regulated  (Eigure  5D).  Several  genes  were 
validated  by  quantitative  PCR  in  a  new  set  of  biological 
replicates  of  468-ERj8#18  cells  treated  with  vehicle  or 


Dox  followed  by  DMSO  or  E2  for  4  hours.  Two  ligand- 
independent  genes  were  confirmed  ( CTGF  and  C3;  Eigure 
6A)  and  three  ligand-dependent  genes  were  confirmed 
{SDCl,  CDHl,  and  DKKl;  Eigure  6B).  Eour  genes  regu¬ 
lated  in  both  a  ligand-dependent  and  a  ligand-indepen¬ 
dent  manner  were  also  validated  (WNT4,  THRSP, 
CABLESl,  and  S100A7;  Eigure  6C).  Importantly,  the  ex¬ 
pression  of  several  ERjS  target  genes  was  assessed  in 
MDA-MB-468-Tet  cells  treated  with  vehicle  or  Dox,  and 
none  of  the  genes  were  regulated  by  Dox  in  the  absence  of 
ERjS  expression  (Supplemental  Eigure  5).  In  addition,  the 
ERjS-selective  ligand  ERB-041  was  found  to  regulate  ERjS 
target  genes  similar  to  E2  (Eigure  7,  A  and  B). 

ER/S  target  genes  are  enriched  in  pathways 
regulating  cell  proliferation,  differentiation, 
and  cell  cycle 

After  globally  identifying  the  ERjS  target  genes  in  ERa- 
negative  TNBC  cells,  we  conducted  a  pathway  analysis  of 
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Table  1.  Most  Highly  Up-regulated  and 
Down-regulated  Genes  in  the  DMSO  vs 

DMS0-l-ERj3  Comparison 

Up-regulated  Genes 

Down-regulated  Genes 

Gene 

Fold 

Fold 

Symbol 

Change 

Gene  Symbol 

Change 

ESR2 

92.93 

TMEM189- 

0.25 

UBE2V1 

FGF19 

44.42 

ACERl 

0.26 

PDZKl 

13.76 

NOV 

0.32 

SCGB2A2 

10.70 

ADAMTS5 

0.41 

THRSP 

9.25 

CTGF 

0.43 

FPHA8 

8.11 

NDRG1 

0.49 

RBP3 

7.55 

KDR 

0.52 

CA12 

6.28 

KRT6B 

0.54 

CYP2A13 

5.60 

NTRK2 

0.55 

MAPK4 

5.04 

KRT14 

0.57 

CD34 

4.77 

IGFBP3 

0.57 

KCNK15 

4.73 

CLIC5 

0.58 

CYP2B7P1 

4.10 

MAOA 

0.59 

C3 

3.69 

FOXOl 

0.59 

ST8SIA6 

3.57 

SYTL4 

0.59 

S100A7A 

3.49 

KRT16 

0.61 

MGAT3 

3.42 

GASP  14 

0.61 

S100A7 

3.39 

PPP1R3C 

0.62 

ANPFP 

3.31 

SMOG! 

0.63 

TMFM229B 

3.04 

KRT17 

0.63 

the  total  ERjS  target  genes  identified  in  comparison  IV 
(Figure  5A)  using  IPA.  These  ERjS  target  genes  were  en¬ 
riched  in  several  molecular  and  cellular  functions,  includ¬ 
ing  cell  movement,  cell  death  and  necrosis,  and  cell  dif- 

Table  2.  Most  Highly  Up-regulated  and 
Down-regulated  Genes  in  the  DMSO-EERjS  vs 

E2-l-ER^  Comparison 

Up-regulated  Genes 

Down-regulated  Genes 

Fold 

Gene 

Fold 

Gene  Symbol 

Change 

Symbol 

Change 

FGR3 

43.43 

ADAMTS5 

0.13 

OTOF 

34.18 

CHRMl 

0.25 

SHISA2 

26.24 

ADAMTSl 

0.27 

ASB2 

18.22 

LOCI  52225 

0.29 

PTH1R 

17.93 

STONl 

0.33 

CYTH4 

15.16 

SOX11 

0.34 

LOCI  00507584 

14.34 

RASSFIO 

0.36 

BSND 

14.19 

KLHL38 

0.39 

TMIF 

10.74 

ARL4A 

0.40 

ENPP2 

10.51 

ARID5B 

0.41 

L0XL4 

9.85 

ABCAl 

0.42 

IP6K3 

8.89 

JAGl 

0.42 

PDK4 

8.80 

VAV3 

0.42 

FAM25A 

8.65 

FAM84A 

0.44 

HAND! 

8.13 

SERTAD4 

0.44 

A4GALT 

8.13 

PPP1R3C 

0.44 

WISP2 

7.88 

RDHIO 

0.44 

VWF 

7.67 

PPP1R3C 

0.62 

THRSP 

7.65 

SMOCl 

0.63 

CABP7 

6.75 

KRT17 

0.63 

Table  3.  Most  Highly  Up-regulated  and  Down- 
regulated  Genes  in  the  E2  vs  E2-I-ER/3  Comparison 


Up-regulated  Genes 

Down-regulated  Genes 

Gene 

Symbol 

Fold 

Change 

Gene 

Symbol 

Fold 

Change 

IS  Ml 

99.57 

ADAMTS5 

0.06 

FGF19 

95.53 

STONl 

0.15 

EGR3 

85.29 

DKKl 

0.16 

WISP2 

69.21 

ADAMTSl 

0.16 

ESR2 

65.01 

CHRMl 

0.20 

THRSP 

54.88 

HRCTl 

0.23 

CYTH4 

53.59 

GRIN2A 

0.23 

OTOF 

51.63 

KLHL38 

0.24 

ENPP2 

39.12 

NOV 

0.24 

SCGB3A1 

36.32 

PPP1R3C 

0.25 

CD34 

29.79 

SOXll 

0.25 

PDLIM4 

nil 

IL8 

0.26 

ANGPTL2 

26.90 

FGFl 

0.26 

PGLYRP2 

26.83 

RASSFIO 

0.27 

LOC727710 

25.48 

FOXQl 

0.28 

SHISA2 

25.25 

SHROOM2 

0.28 

CRB2 

25.05 

ARL4A 

0.28 

TMIF 

24.75 

ACERl 

0.28 

FOXNl 

24.04 

F AM  105 A 

0.29 

MGAT3 

23.58 

IGFBP3 

0.29 

ferentiation  (Table  4).  In  addition,  genes  involved  in  cell 
proliferation  and  cell  morphology  were  enriched  in  the 
ER/3  target  genes.  Similar  molecular  and  cellular  func¬ 
tions  were  enriched  in  genes  identified  in  comparisons  II 
and  III  (Figure  5A),  representing  potential  ligand-inde¬ 
pendent  and  ligand-dependent  targets,  respectively  (Sup¬ 
plemental  Materials  and  Methods).  The  top  pathways 
enriched  in  the  total  ERjS  target  gene  set  included  the 
Wnt//3-catenin  pathway  and  the  Gl/S  cell  cycle  check¬ 
point  regulation  pathway  (Table  5).  Indeed,  several  genes 
involved  in  Wnt  signaling  were  verified  as  ERjS  target 
genes,  including  DKKl,  WNT4,  and  CDHl  (Figure  6,  B 
and  C).  In  addition,  the  gene  encoding  the  cyclin-depen- 
dent  kinase  inhibitor  p21,  CDKNIA,  was  found  to  be 
up-regulated  by  ERjS  expression  and  E2  treatment  (Figure 
7A).  The  pathway  analysis  results  support  the  potential 
for  ERjS  to  regulate  key  cellular  process  that  may  ulti¬ 
mately  play  a  role  in  growth  inhibition. 

ER/S  target  genes  are  regulated  in  additional 
breast  cancer  cell  lines  and  correlate  with  ER/S 
expression  in  TNBCs 

To  confirm  that  a  common  set  of  ERjS  target  genes  are 
shared  among  TNBC  cells,  the  genes  identified  in  MDA- 
MB-468-ERjS  cells  were  compared  to  those  identified  in 
other  TNBC  cell  lines.  First,  the  gene  list  was  compared  to 
that  obtained  by  Secreto  and  coworkers  (15)  in  which 
Hs578T  breast  cancer  cells  with  inducible  ER/S  expres¬ 
sion  were  used  to  identify  ligand-dependent  ER/S  target 
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Figure  6.  Validation  of  ER/3  target  genes  in  468-ERj3#18  cells.  468- 
ER(3#18  cells  were  treated  with  vehicle  or  Dox  for  48  hours  followed 
by  treatment  with  DMSO  or  10  nM  E2  for  4  hours.  Quantitative  PCR 
was  performed  as  described  in  the  Materials  and  Methods.  (A) 
Expression  of  ESR2  and  2  ligand  independent  target  genes:  CTGF  and 
C3.  (B)  Expression  of  three  ligand-dependent  target  genes:  SDC1 , 
CDHl ,  and  DKK1 .  (C)  Expression  of  four  genes  regulated  in  both 
ligand-independent  and  ligand-dependent  manners:  WNT4,  THRSP, 
CABLESl ,  and  S100A7.  Data  represent  the  mean  and  SEM  of  three 
independent  replicates.  *,  P  <  .05  compared  to  DMSO  control; 

#,  P  <  .05  compared  to  DMSO  +  Dox  control. 


genes  after  24  hours  of  E2  treatment.  Although  the  au¬ 
thors  used  a  targeted  microarray  to  assess  changes  in  gene 
expression,  27  genes  were  found  to  be  commonly  regu¬ 
lated  in  both  MDA-MB-468-ER/3  and  l-Is578T-ERj8  cells 
(Table  6).  Most  of  these  genes  were  regulated  in  the  same 
direction,  with  the  exceptions  of  BMPS,  CEBPD,  and 
PMAIPI.  Eour  ERj8  target  genes  were  then  verified  in 
another  TNBC  cell  line,  EICC1143,  infected  with  retro¬ 
virus  encoding  green  fluorescent  protein  (GEP)  as  a  con¬ 
trol  or  full-length  ERjS  followed  by  treatment  with 
DMSO  or  E2  for  4  hours  (Eigure  7C).  Indeed,  all  four 
genes  including  the  cell  cycle  regulator  CDKNIA  were 


A 

■  DMSO  S  E2  □  DMSO+Dox  Q  E2+Dox 


B 


C 


GFP/DMSO  S  GFP/E2  □  ERp/DMSO  E  ERp/E2 


Figure  7.  ER^  expression  regulates  common  target  genes  in 
HCC1 143  TNBC  cells.  (A)  Quantitative  PCR  validation  of  SEMA3B, 
MXD1 ,  and  CDKNIA  in  468-ERj3#18  cells  treated  with  vehicle  or  Dox 
for  24  hours  followed  by  treatment  with  DMSQ  or  10  nM  E2  for  4 
hours.  Data  represent  the  mean  and  SEM  of  three  independent 
replicates.  *,  P  <  .05  compared  to  DMSQ  control;  #,  P  <  .05 
compared  to  DMSO  +  Dox  control.  (B)  Quantitative  PCR  validation  of 
ESR2,  SEMA3B,  MXD1 ,  and  CDKNIA  in  HCC1 143  TNBC  cells  infected 
with  retrovirus  for  GEP  or  ERj3  expression  followed  by  treatment  with 
DMSO  or  1 0  nM  E2  for  4  hours.  The  data  represent  the  mean  and  SEM 
of  three  replicates.  *,  P  <  .05  compared  to  GFP/DMSO  control;  #,  P  < 
.05  compared  to  ERjS/DMSO  control.  (C)  Quantitative  PCR  validation  of 
SEMA3B,  MXD1,  and  CDKNIA  in  468-ERj3#18  cells  treated  with 
vehicle  or  Dox  for  24  hours  followed  by  treatment  with  DMSO  or  100 
nM  ERB-041  for  4  hours.  Data  represent  the  mean  and  SEM  of  three 
independent  replicates.  *,  P  <  .05  compared  to  DMSO  control;  #,  P  < 
.05  compared  to  DMSO  +  Dox  control. 


regulated  in  a  similar  manner  in  EICC1143  when  ERjS 
was  expressed  and  activated  by  E2  treatment.  These  data 
suggest  that  ERjS  regulates  multiple  common  targets 
across  TNBC  cell  lines. 

The  ERjS  target  genes  identified  in  MDA-MB-468-ERjS 
cells  were  next  compared  to  ERjS  target  genes  identified 
by  Grober  and  colleagues  (30)  in  ERa-positive  MCE7 
cells  engineered  to  express  full-length  ERjS.  Using  chro¬ 
matin  precipitation  followed  by  deep  sequencing  for  ERjS 
binding  sites  and  gene  expression  microarrays,  they  found 
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Table  4.  Top  Biological  Functions  in  the  E2  vs  E2  +  ERj3 
Gene  Set 


Top  Biological  Functions, 
Molecular  and  Cellular 
Functions 

P  Value 

No.  of 
Molecules 

Cell  movement 

Cell  movement 

1.21E-18 

166 

Migration  of  celis 

2.73E-18 

154 

Ceil  death  and  survival 

Necrosis 

2.19E-17 

201 

Celi  death 

2.66E-16 

236 

Ceilular  development 

Differentiation  of  celis 

6.66E-17 

129 

Proiiferatlon  of  tumor  cell  lines 

1.54E-12 

143 

Ceilular  growth  and  proliferation 

Proliferation  of  tumor  cell  lines 

1.54E-12 

143 

Proliferation  of  cells 

1.70E-12 

249 

Cellular  morphology 

Formation  of  lamellipodia 

4.56E-09 

22 

Collapse  of  growth  cone 

2.31  E-07 

12 

424  “primary”  ERjS  target  genes  that  were  differentially 
expressed  when  ER/3  was  expressed  and  contained  ERj8 
binding  sites  within  10  kb  of  the  gene.  Of  the  930  ER/3 
target  genes  identified  in  MDA-MB-468-ER/3  cells  (com¬ 
parison  TV,  Eigure  5A),  99  genes  (10.6%)  were  also  iden¬ 
tified  as  primary  ER/3  targets.  Eifty-two  of  these  genes 
were  differentially  regulated  in  cells  expressing  ER/3  and 
not  in  parental  MCE7  cells  (Supplemental  Materials  and 
Methods).  These  results  indicate  that  despite  the  differ¬ 
ences  in  cell  lines  (ie,  ERa-negative  basal-like  cells  versus 
ERa-positive  luminal  cells),  ER/3  can  regulate  common 
target  genes.  Most  of  these  genes  were  also  regulated  in 
the  absence  of  ERj8  or  contained  ERa  binding  sites,  but 
three  genes  were  found  to  contain  only  ERjS  binding  sites 
and  were  only  regulated  when  ERjS  was  expressed: 
CCDC103,  IGSF9,  and  P2RY6.  Although  CCDC103 
was  only  identified  in  comparison  IV  (Eigure  5A),  IGSF9 
was  found  to  be  a  ligand-dependent  target  (comparison 
III)  and  P2RY6  was  a  ligand-dependent  target  (compari¬ 
son  II)  in  MDA-MB-468-ERjS  cells. 


Table  5.  Top  Pathways  in  the  E2  vs  E2-FERjS  Gene  Set 


Top  Canonical  Pathways 

P  Value 

Ratio 

Flepatic  fibrosis/hepatic  stellate  cell 
activation 

7.24E-07 

0.16 

Axonal  guidance  signaling 

2.57E-05 

0.09 

Wnt/j3-catenin  signaling 

1.9  IE-04 

0.12 

Role  of  IL-17A  in  psoriasis 

3.63E-04 

0.39 

Aryl  hydrocarbon  receptor  signaling 

6.17E-04 

0.11 

Molecular  mechanisms  of  cancer 

8.51  E-04 

0.09 

Glioblastoma  multiforme  signaling 

1 .66E-03 

0.10 

Fluman  embryonic  stem  cell  pluripotency 

1 .82E-03 

0.10 

Coagulation  system 

2.04E-03 

0.18 

Glutamate  receptor  signaling 

2.45E-03 

0.13 

Cell  cycle;  G1/S  checkpoint  regulation 

3.55E-03 

0.14 

To  validate  the  ERjS  target  genes  in  TNBCs,  a  correla¬ 
tion  analysis  was  performed  using  gene  expression  data 
from  43  TNBCs  available  through  TCGA  (26).  Of  the 
930  ERjS  target  genes  identified  in  MDA-MB-468-ERjS 
cells  (comparison  TV,  Eigure  5A),  109  genes  were  signifi¬ 
cantly  correlated  with  ESR2  expression  {q  <  .05,  denoted 
by  points  in  Eigure  8;  a  complete  list  is  provided  in  Sup¬ 
plemental  Eile  2).  Eifteen  of  these  were  found  to  be  ligand- 
independent  in  MDA-MB-468-ERjS  (identified  compari¬ 
son  II,  Eigure  5A),  and  58  of  them  were  regulated  in  a 
ligand-dependent  manner  (identified  in  comparison  III, 
Eigure  5A).  Six  of  these  genes  were  also  identified  in 
Hs578T-ERj8  cells  (highlighted  in  red  in  Eigure  8), 
thereby  validating  that  these  genes  are  likely  ERjS  target 
genes  in  TNBCs.  However,  two  genes  (ADAMTSl  and 
ADAMTS5)  showed  a  positive  correlation  with  ERjS  ex¬ 
pression  but  were  down-regulated  in  both  cell  lines,  sug¬ 
gesting  that  ERjS  target  gene  regulation  in  breast  cancer 
samples  is  not  completely  recapitulated  by  the  cell  culture 
models.  In  addition,  two  of  the  genes  identified  as  primary 
ERjS  target  genes  by  Grober  et  al  (30)  were  significantly 
correlated  with  ESR2  expression  in  the  TCGA  cohort 
(highlighted  in  green  in  Eigure  8),  although  the  correla¬ 
tion  coefficients  were  quite  low.  This  could  be  due  to  the 
difference  in  breast  cancer  subtypes:  Grober  and  col¬ 
leagues  used  an  ERa-positive  cell  line  to  identify  ERjS 
target  genes,  while  only  TNBCs  were  included  in  the 
TCGA  analysis.  Overall,  these  results  validate  that  a  set  of 
ERjS  target  genes  identified  in  cell  culture  models  are  cor¬ 
related  with  ERjS  expression  in  TNBCs. 

Discussion 

Because  ERjS  is  expressed  in  50%  to  90%  of  ERa-nega- 
tive  breast  cancers  (1-4)  and  its  expression  has  been  as¬ 
sociated  with  a  better  prognosis  in  TNBCs  (2),  ERjS  may 
be  a  therapeutic  target  for  a  subset  of  TNBCs.  The  results 
of  this  study  support  this  hypothesis  and  confirm  the 
growth  inhibitory  effects  of  ERjS  expression  and  activa¬ 
tion  in  TNBC  cells.  However,  the  requirement  of  ligand  to 
stimulate  ERjS’s  effects  is  likely  cell  context  specific.  Pre¬ 
viously,  ERjS  expression  was  found  to  inhibit  ERa-nega- 
tive  breast  cancer  cell  growth  independent  of  ligand  in 
MDA-MB-231  and  SKBR3  cells  (12,  13).  In  Hs578T 
TNBC  cells,  ERjS  expression  inhibited  cell  growth  in  a 
ligand-dependent  manner  through  a  G1  cell  cycle  arrest 
(15).  In  culture,  ERjS  expression  also  induced  a  G1  arrest 
in  MDA-MB-468  cells,  but  this  effect  was  both  E2-in de¬ 
pendent  and  -dependent.  One  way  that  ERjS  may  regulate 
the  cell  cycle  is  through  up-regulation  of  the  cyclin-depen- 
dent  kinase  inhibitor  p21  (encoded  by  CDKNIA),  which 
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Table  6.  Common  ERjS  Target  Genes  Identified  in  Hs578T-ERjS  (1 5)  and  468-ERjS#18 

Gene 

Symbol 

Gene  Name 

Fold  Change 
MDA468-ER/3 

Fold  Change 
HS578T-ER/3 

ADAMTSl 

ADAM  metallopeptidase  with  thrombospondin  type  1  motif,  1 

0.2 

0.4 

ADAMTS5 

ADAM  metallopeptidase  with  thrombospondin  type  1  motif,  5 

0.1 

0.4 

BMP5 

Bone  morphogenetic  protein  5 

0.5 

2.6 

C3 

Complement  component  3 

4.4 

12.8 

CAT  2 

Carbonic  anhydrase  XII 

10.5 

2.9 

CD34 

CD34  molecule 

29.8 

2.1 

CEBPD 

CCAAT/enhancer  binding  protein  (C/EBP),  8 

0.5 

2.7 

CTGF 

Connective  tissue  growth  factor 

0.3 

0.4 

DKKl 

Dickkopf  homolog  1  {Xenopus  laevis) 

0.2 

0.3 

EHD2 

EH-domain  containing  2 

2.9 

2.3 

ENPP2 

Ectonucleotide  pyrophosphatase/phosphodiesterase  2  (autotaxin) 

39.1 

2.0 

EUT8 

EucosyltransferaseS  (a  [1,6]  fucosyltransferase) 

2.1 

2.1 

IGEBP3 

Insulin-like  growth  factor  binding  protein  3 

0.3 

0.4 

IGEBP5 

Insulin-like  growth  factor  binding  protein  5 

4.6 

2.5 

JAG1 

Jagged  1  (Alagille  syndrome) 

0.3 

0.1 

LM02 

LIM  domain  only2  (rhombotin-like  1) 

1.5 

3.3 

PDE4B 

Phosphodiesterase  4B,  cAMP-specific 

2.0 

2.3 

PFKEB3 

6-phosphofructo-2-kinase/fructose-2,6-biphosphatase  3 

4.5 

2.1 

PI  Ml 

Pim-1  oncogene 

2.2 

2.3 

PLACl 

Placenta-specific  1 

2.2 

3.1 

PLAU 

Plasminogen  activator,  urokinase 

0.5 

0.4 

PMAIPl 

Phorbol-12-myristate-1 3-acetate-induced  protein  1 

2.1 

0.5 

SEMA3B 

Serna  domain,  immunoglobulin  domain  (Ig),  (semaphorin)  3B 

7.2 

3.6 

SFRPl 

Secreted  frizzled-related  protein  1 

1.8 

2.1 

TH 

Tyrosine  hydroxylase 

17.7 

2.2 

TPD52L1 

Tumor  protein  D52-like  1 

1.6 

2.0 

WISP2 

WNT1  inducible  signaling  pathway  protein  2 

69.2 

3.0 

regulates  the  progression  from  G1  to  S  phase.  Indeed,  one 
of  the  top  pathways  identified  in  ERjS  target  genes  was  the 
Gl/S  checkpoint  regulation  pathway  (Table  5),  which  in¬ 
cludes  CDKNIA.  In  MDA-MB-231  TNBC  cells,  ERjS 
increased  the  expression  of  p21  via  interaction  with  an¬ 
other  transcription  factor,  inhibitor  of  differentiation- 1 
(14).  Euture  studies  are  warranted  to  determine  if  inhib¬ 
itor  of  differentiation- 1  plays  a  role  in  the  regulation  of 
CDKNIA  in  MDA-MB-468  cells. 

ERjS-mediated  growth  inhibition  may  also  depend  on  the 
cellular  environment.  ERj8  was  found  to  elicit  growth  inhib¬ 
itory  effects  in  vitro  and  in  vivo  in  MDA-MB-468-ERj8  cells. 
In  contrast,  in  mammary  epithelial  and  breast  cancer  cells 
that  endogenously  coexpress  ERa  and  ERjS,  ERjS-selective 
ligands  were  found  to  induce  apoptosis  in  vitro  while  pro¬ 
moting  tumor  growth  in  vivo  (31)  by  coincubation  with  a 
basement  membrane  extract,  suggesting  that  the  cellular  en¬ 
vironment  can  modulate  ERjS  activity.  Inhibition  of  extra¬ 
cellular  signal-regulated  kinase  1  and  2  (ERKl/2)  could  re¬ 
store  ERjS  growth  inhibitory  effects  (31).  It  appears  that 
cellular  context  modulates  ERjS  growth  inhibitory  effects. 
Either  the  coexpression  of  ERa  or  the  activation  of  kinase 
pathways  in  the  mammary  epithelial  cell  model  may  account 
for  the  differential  in  vivo  effects  of  ERjS  in  the  MDA-MB- 
468  orthotropic  xenograft  model. 


Like  ERa,  ERjS  has  two  activation  functions  (AE-1  and 
AE-2)  that  confer  transcriptional  activity.  The  AE-1  can 
mediate  transcription  independent  of  ligand,  and  the 
AE-2  functions  in  a  ligand-dependent  manner  (32).  Al¬ 
though  the  AE-1  of  ERjS  was  found  to  have  negligible 
activity  in  a  Gal4  assay  (33),  ERjS  was  found  to  interact 
with  DNA  in  a  ligand-independent  manner  (34).  In  addi¬ 
tion,  MAP  kinase  has  been  shown  to  phosphorylate  the 
N-terminal  activation  function  (AE-1)  of  ERjS,  thereby 
promoting  the  recruitment  of  steroid  receptor  coactivator 
1,  which  confers  ligand-independent  transcriptional  ac¬ 
tivity  (35).  In  MDA-MB-468  cells,  ERjS  appears  to  have 
both  ligand-independent  and  -dependent  activity  in  re¬ 
gard  to  target  gene  expression.  Interestingly,  most  target 
genes  were  regulated  in  response  to  E2  treatment,  suggest¬ 
ing  strong  ligand-dependent  activity  in  these  cells.  In  con¬ 
trast,  most  ERjS  target  genes  were  regulated  independent 
of  ligand  in  U20S  osteosarcoma  cells  (36).  Therefore,  the 
ligand-independent  activity  of  ERjS  is  likely  the  result  of 
the  availability  of  cofactors,  the  activation  of  signaling 
pathways,  and  the  promoter  context. 

Both  the  ligand-dependent  and  -independent  ERjS  tar¬ 
get  genes  identified  in  MDA-468-ERjS  cells  were  enriched 
in  genes  that  regulate  cell  death,  proliferation,  movement, 
and  morphology.  ERjS  target  genes  identified  after  24 
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Figure  8.  ERjS  target  genes  are  correlated  with  ERj3  expression  in  a 
cohort  of  TNBCs.  Gene  expression  data  from  the  TCGA  (26)  was  used 
to  perform  a  correlation  analysis  of  ERj3  target  genes  and  ERj3 
expression.  Of  the  930  genes  identified  468-ERj3#1 8  cells  (comparison 
IV,  Figure  5A),  109  genes  were  significantly  correlated  with  ER/3 
expression  (q  <  .05).  The  correlation  coefficient  versus  the  fold  change 
observed  in  the  RNA-seq  results  for  each  gene  is  shown  in  the 
scatterplot.  The  gene  symbols  of  those  genes  with  the  highest 
correlation  coefficients  are  shown  on  the  plot.  Genes  in  red  are  those 
that  were  also  identified  in  Hs578T-ERj3  cells  (1 5),  and  genes  in  green 
are  those  that  were  identified  in  l\/ICF7-ER/3  cells  (30). 

hours  of  E2  treatment  in  HEK293  cells  were  found  to  be 
enriched  in  many  of  the  same  pathways  (37).  A  pathway 
analysis  of  total  ERjS  target  genes  revealed  significant  en¬ 
richment  of  several  pathways,  including  the  Wnt/jS- 
catenin  pathway.  Of  the  genes  involved  in  the  Wnt/jS- 
catenin  pathway,  ERjS  expression  and  E2  treatment 
resulted  in  the  up-regulation  of  WNT4  (a  noncanonical 
Wnt  ligand)  and  CDHl  (a  jS-catenin  interacting  protein), 
as  well  as  the  down-regulation  of  DKKl  (an  inhibitor  of 
Wnt  signaling).  Down-regulation  of  DKKl  by  ERjS  ex¬ 
pression  and  estrogen  treatment  was  also  observed  in 
l-Is578T  breast  cancer  cells  (15)  and  293T  kidney  cells 
(37).  Collectively,  the  gene  expression  data  suggest  that 
ERjS  expression  may  affect  the  Wnt  signaling  pathway  in 
breast  cancer  cells.  Although  activation  of  Wnt  signaling 
can  transform  human  mammary  epithelial  cells  (38),  a 
recent  study  has  shown  that  paracrine  Wnt  signaling  in 
fibroblasts  can  either  promote  or  inhibit  breast  cancer  cell 
growth  in  mice  (39).  ERjS  has  been  detected  in  both  mam¬ 
mary  epithelial  and  stromal  cells  (40)  and  may  therefore 
affect  Wnt  signaling  in  a  cell  type-  and  cell  context-depen¬ 
dent  manner. 

Given  the  heterogeneous  nature  of  TNBCs  (41),  iden¬ 
tification  of  common  ERjS  target  genes  across  all  TNBCs 
would  be  challenging.  Recently  TNBC  cell  lines  were  clas¬ 
sified  into  basal-like,  mesenchymal-like,  and  luminal-an¬ 
drogen  receptor  subtypes  by  correlation  with  expression 


profiles  of  large  breast  cancer  data  sets  (42).  Although  the 
RNA-seq  was  performed  in  basal-like  MDA-MB-468 
cells,  we  verified  ERjS  target  gene  expression  in  EIs578T, 
a  mesenchymal-like  cell  line,  and  EICC1143,  another  bas- 
al-like  cell  line.  Despite  differences  in  the  cell  line  and 
experimental  design,  27  genes  identified  in  MDA-MB- 
468-ERj3  cells  were  found  to  be  ERjS  target  genes  in 
EIs5  78T-ERj8  cells  treated  with  E2  for  24  hours  (15).  Eour 
ERjS  target  genes  selected  for  verification  showed  the  sim¬ 
ilar  regulation  by  ERjS  and  E2  in  EICC1143  and  MDA- 
MB-468  cells  (Eigure  7C).  In  addition,  an  ERjS  target  gene 
identified  in  MDA-MB-468-ERjS  cells,  S100A7,  was 
found  to  be  an  ERjS-specific  target  gene  in  MCE7  ERa- 
positive  breast  cancer  cells  with  inducible  expression  of 
ERjS  and  was  associated  with  ERjS  in  ERa-negative  breast 
tumors  (43).  Several  ERjS  targets  were  also  found  to  be 
primary  ERjS  target  genes  in  MCE7  cells  (30)  and  are 
likely  direct  targets  of  the  receptor  in  MDA-MB-468-ERjS 
cells.  Most  importantly,  109  ERjS  target  genes  were  cor¬ 
related  with  ERjS  expression  in  a  cohort  of  43  TNBCs, 
including  six  genes  that  were  also  identified  in  Hs578T 
TNBC  cells.  Thus,  the  genes  identified  in  this  study  pro¬ 
vide  the  foundation  for  identifying  common  target  genes 
of  ERjS  functionality  across  ERa-negative  breast  cancers, 
including  TNBCs. 

Overall,  this  study  confirms  the  growth  inhibitory  ef¬ 
fects  of  ERjS  in  TNBC  cells,  both  in  vitro  and  in  vivo,  and 
provides  the  most  comprehensive  identification  of  ERjS 
target  genes  in  TNBC  to  date.  These  data  also  suggest  that 
ERjS  selective  ligands  may  be  useful  for  targeting  ERjS  in 
a  subset  of  ERa-negative  breast  cancers,  including 
TNBCs.  ERjS-selective  ligands  could  promote  the  inhibi¬ 
tory  effects  of  ERjS  while  avoiding  the  proliferative  effects 
mediated  by  ERa.  ERjS  selective  ligands,  such  as  ERB- 
041,  have  been  found  to  be  well-tolerated  (44),  support¬ 
ing  the  use  of  these  compounds  clinically.  In  addition,  the 
identification  of  ERjS  target  genes  provides  a  foundation 
for  identifying  biomarkers  of  ERjS  functionality  in 
TNBCs.  Collectively,  these  data  support  the  hypothesis 
that  ERjS  may  be  a  therapeutic  target  in  a  subset  of 
TNBCs  and  suggest  that  ERjS  may  be  targeted  with  ERjS- 
selective  ligands  that  can  stimulate  its  growth  inhibitory 
effects. 
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